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SPHEROIDAL WAVE OPERATORS WITH COMPLEX 
ASPHERICAL PARAMETER 
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Abstract. A family of spectral decompositions of the spin-weighted spheroidal 
wave operator is constructed for complex aspherical parameters with bounded imag¬ 
inary part. As the operator is not symmetric, its spectrum is complex and Jordan 
chains may appear. We prove uniform upper bounds for the length of the Jor¬ 
dan chains and the norms of the idempotent operators mapping onto the invariant 
subspaces. The completeness of the spectral decomposition is proven. 
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1. Introduction and Statement of Results 


The spin-weighted spheroidal wave equation arises in the study of electromagnetic, 
gravitational and neutrino-field perturbations of rotating black holes when separating 
variables in the so-called Teukolsky master equation (see la |22] or the survey pa¬ 
per m)- In the spin-weighted spheroidal wave equation, the spin of the wave enters as 
a parameter s G {0, 1, |, 2,...}. We are mainly interested in the cases s = 1 of an 

electromagnetic and s = 2 of a gravitational field. If s is an integer, the spin-weighted 
wave equation is the eigenvalue equation 


A'i’ = X'i’, 


( 1 . 1 ) 


where the spin-weighted spheroidal wave operator A is an elliptic operator with smooth 
coefficients on the unit sphere More specifically, choosing polar coordinates i9 G 
(0, tt) and (p G [0,27r), we have (see for example [21]) 


. d . 2 a 9 

A = — — -^ sm '& 


d cos'd 


d cos'd 


H-\—fu sin^ t?-|-- scos'd') . 

sin^ 7? V dip / 


Here D G C is the aspherical parameter. In the special case D = 0, we obtain the spin- 
weighted Laplacian on the sphere, whose eigenvalues and eigenfunctions can be given 
explicitly m- In the case s = 0 and D / 0, one gets the spheroidal wave operator 
mu)- Setting D = 0 and s = 0, one simply obtains the Laplacian on the sphere. 
We consider A on the Hilbert space !K = with domain of definition Ti[A) = 

C°°{S‘^). We remark that A clearly is an elliptic operator on the sphere. However, 
even in the case s = 0 and for real D, in general there is no Riemannian metric on the 
sphere which realizes the spheroidal wave operator as the Laplace-Beltrami operator. 
Thus the spheroidal wave operator cannot be be identified with the Laplace-Beltrami 
operator on a Riemannian manifold. Eor general spin, this means in particular that 
the methods used for spin-weighted spherical harmonics in m Section 4.15] do not 
seem to generalize to the spheroidal situation. 
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As the spin-weighted spheroidal wave operator is axisymmetric, we can separate out 
the (/9-dependence with a plane wave ansatz, 

ip) = 0(i9) with A: G Z . 


Then A becomes the ordinary differential operator 


Ak := -- 


d 


sin^ 'd ■ 


d 


1 


-|- — (n sim d + k — s cos i?)" 


( 1 . 2 ) 


d cos I? d cos 'd ' sin^ i? 

This operator acts on the vectors in !K with the prescribed (/)-dependence, which we 
denote by "Kk, 

Ak := L‘^{S^) n 0(i9) | 0 : (0, vr) ^ C} . 

The domain of definition reduces to 


V{Ak) = C°^{S^)nAk. 

The Hilbert space Ak can be identified with 

Ak = L^((0,7r),sini?(ii?) . 

Also, one can consider Ak as an ordinary differential operator on this Hilbert space, 
for example with the domain of definition (^““((0,7r))nL^((0, tt). However, when doing 
so, one still needs to specify boundary conditions at i? = 0, vr. As will be explained in 
detail in Section [2] below, the correct boundary conditions are that the limits 

lim 0(i?) must exist . (1.3) 

i9^0,7r 

In this formulation as a pure ODE problem, the spheroidal wave equation ()1.2I) can 
also be used in the case of half-integer spin (to describe neutrino or Rarita-Schwinger 
fields), if k is chosen to be a half-integer. Thus in what follows, we fix the parameters s 
and k such that 

2s G No and A; — s G Z. (1.4) 

We are interested in the case that D is complex. Then the potential in (jl.2ll is 
complex, so that the operator Ak is not a symmetric operator on Ak- As a consequence, 
the spectral theorem in Hilbert spaces does not apply. The spectrum will in general be 
complex. Moreover, the operator need not be diagonalizable, because Jordan chains 
may form. The main task of the present paper is to control the length of these Jordan 
chains to obtain a decomposition of the Hilbert space into invariant subspaces of Ak- 
This is our main result: 


Theorem 1.1. For any s and k in the range dm) and any c > 0, we let U G C be 
the strip 

|ImD| < c . (1.5) 

Then there is a positive integer N and a family of bounded linear operators Qn{'^) 
on Ak defined for all n gNU {0} and Tl gU with the following properties: 

(i) The image of the operator Qq is an N-dimensional invariant subspace of Ak- 

(ii) For every n > 1, the image of the operator Qn is an at most two-dimensional 
invariant subspace of Ak- 

(hi) The Qn are uniformly bounded in TifAk), i.e. for all n G N U {0} and Q G U, 

WQnW < C2 ( 1 . 6 ) 

for a suitable constant c^ = C 2 (s, fe, c) (here || • || denotes the sup-norm on Ak). 
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(iv) The Qn are idempotent and mutually orthogonal in the sense that 

Qn Qn' ^n,n' Qn far all 71, 71 G U {0} . 

(v) The Qn are complete in the sense that for every Q ^ U, 

oo 

J2Qn = l ( 1 - 7 ) 

n =0 

with strong convergence of the series. 

Note that the operators Qn are in general not symmetric (i.e. Q* 7 ^ Qn)- This corre¬ 
sponds to the fact that for non-symmetric operators, the eigenvectors corresponding 
to different eigenvalues are in general not orthogonal. 


2. Reformulation as a Sturm-Liouville Problem 


We first bring the operator (|1.2p to the standard Sturm-Liouville form (for more 
details see m Section 2]). To this end, we first write the operator in the variable u = 
id G (0,7r), 

A i d , d 1 / .2 7 \ 2 

Ak =-sinu -—|- 7 ^— (if sm u + k — s cos u) . 

sm u du du sin u 

Introducing the function (j) by 

(j) = Vsin u 0 , 

we get the eigenvalue equation 

Hf) = X(f> , 

where H has the form of a one-dimensional Hamiltonian 

d^ 




where W is the complex potential 


1 cos^ 71 1 1 /r-. 2 7 

Vv = — T —T-- -|———(iZ sm u + k — s cos uj" 


4 sin^ u 2 


sm^ 71 


= sin^ u + [ k'^ + s‘^ — 


1 


1 


sm 71 


-h 2nk -s^ - 


— 2sR cos 71 — 2sk 


cos 71 


sin^ 71 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 


For what follows, it is usually most convenient to write (12.2p as the the Sturm-Liouville 
equation 

-^ + 00 = 0, (2.7) 


where V is the potential 

V = sin^ 71-1- (k'^ + s'^ — 


1 


1 


4 / sin^ 71 


— 2sR cos 71 — 2sk 


cos 71 
sin^ 71 




( 2 . 8 ) 


and pL is the constant 

pi = \-29.k +s‘^+ ^. (2.9) 

The transformation (|2.ip from 0 to P becomes a unitary transformation if the integra¬ 
tion measure in the corresponding Hilbert spaces is transformed from sin u du to du. 
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Hence the eigenvalue equation (|l.lll on is equivalent to ()2.7I) on the Hilbert space 
-L^(( 0 , TT),du). 

If the potential V were continuous on the interval [0, vr], we would get a well-defined 
boundary problem by imposing Dirichlet or Neumann or more general mixed boundary 
values at 0 and vr (for details see O Chapter 12]). In our situation, there is the com¬ 
plication that the potential (|2.8p has poles at the boundary points. As a consequence, 
the fundamental solutions will also have singularities, so that it is no longer obvious 
how to introduce suitable boundary conditions. In the case when s is an integer, the 
correct boundary values can be determined by going back to the eigenfunctions on the 
sphere as we now explain. Due to elliptic regularity theory, the eigenfunctions 'k 

of the angular operator are smooth functions on the sphere. Therefore, we ob¬ 
tain (|1.3I1 as a necessary condition. In view of the transformation (|2.1I) . this implies 
that the limits 

lim n “2 cj){u) and lim (vr — u )~2 cj){u) must exist . ( 2 - 10 ) 

U\^0 U/^TT 

These boundary conditions can also be understood by looking at the asymptotics of 
the solutions of (|2.7p near the boundary points. Namely, expanding the potential (|2.8p 
near the boundary points, we obtain 

V ( u ) = ^ ((t - »)" - j) + 0(u-') 

rW = - j) + o(P - ■ 

If the factors k± s are non-zero, the solutions have the asymptotics 

(/>(n) ~ (l-|- 0 (n)) and (/>(n) ~ (vr — u) (l-|- 0 (vr — n)) . ( 2 . 11 ) 

If on the other hand, the factors A: ± s are zero, the asymptotic solutions involve an 
additional logarithm (for detail see [T 2 l Sections 7 and 8 ]), 

(l){u) = Cl y/u + C 2 y/u ^ogu + 0{u) if k = s (2.12) 

4’{u) = Cl y/ir — u -|- C 2 y/n — u log(vr — u) + 0(vr — u) if k = —s . (2-13) 

In each case, the boundary conditions (I2.10h single out one of the two fundamental 
solutions. In this way, the conditions (I2.inh give mathematically reasonable boundary 
conditions. We remark that in the case (|2.1ip . our boundary conditions are equivalent 
to Dirichlet boundary conditions. Alternatively, these boundary conditions could be 
implemented simply by demanding that the eigenfunctions must be square integrable 
(note that, in view of (|1.4p . the parameter k—s is always an integer). In the exceptional 
cases (I2.12P and (|2.13l) . however, both fundamental solutions satisfy Dirichlet boundary 
conditions and are square integrable. Thus in these cases, it is essential to state the 
boundary conditions in the form (I2.10p . 

In order to bring the boundary conditions (I2.10p into a more tractable form, it 
is convenient to work with solutions of the corresponding Riccati equation: For any 
solution (j) of the Sturm-Liouville equation (|2.7I) . the function y := 4>'/4> satisfies the 
corresponding Riccati equation 

y’ = V-y\ (2.14) 

Using the results and methods in uni US], we can construct a solution y of the Riccati 
equation with rigorous error bounds. With this in mind, let us assume that a solution y 
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of the Riccati equation is known. Then a particular solution of the corresponding 
Sturm-Liouville equation is obtained by integration, 

(j){u) = ex.p ( f yY (2.15) 

^ JUQ ^ 

The general solution can be constructed by integrating the equation for the Wronskian. 
Namely, if ^ is another solution of the Sturm-Liouville equation, the Wronskian 

$) := (l)Y - 


is a constant, and thus 


^{u) = 4>{u) 


f <^(^o) _ r 

yHuo) Juo 


w{4>A) 


(2.16) 


4> l { u ) 




In particular, this relation can be used to construct solutions of the Sturm-Liouville 
equation (12.71) which satisfy the boundary conditions ()2.10p . We denote these solutions 
by (j)^ and (where the subscript D refers to “Dirichlet”, and L/R to the left and 
right boundary points at u = 0 and u = tt, respectively). To this end, we let (j)L and (j)R 
be generic solutions which do not satisfy the boundary conditions (I2.10p . i.e. 

y^^-\k-s\ _|_ 0 (ti)) if fc / S 

^/u logw (l -I- 0(ri)) if A: = s 

(tt — tt) 2 “l*'+®l (l -|- 0(7r — u)) if A: 7 ^ —s 

■\/7r — u log(7r — u) (l -|- 0(7r — n)) ii k = —s . 

Then, using (12.161) . the solutions which do satisfy (I2.10p are given (up to irrelevant 
prefactors) by 

4>Liu) = 4 >l{u) [ and 4>^{u) =-4 >r{u) [ (2.17) 

Jo 91 Ju 9r 

If 0 is a solution of (|2.7p subject to the boundary conditions (I2.10p . then this solution 
must be a multiple of both and (JJr. Hence 4Y Yr linearly dependent, and 
their Wronskian vanishes, 

=0. (2.18) 

In this way, we have reformulated the existence problem for solutions satisfying the 
boundary conditions (j2.I0l) in terms of the vanishing of the Wronskian (|2.18l) . More 
generally, the Wronskian can be used to describe the spectrum of the Hamiltonian. 
Namely, we just saw that if the Wronskian vanishes, then there is an eigensolution 
which satisfies the boundary conditions (j2.7p . Conversely, if the Wronskian is non¬ 
zero, we may introduce the Green’s function by 

1 i Yliu) ifu<u' 


s\{u,u ) = 


[(pLiu')4'Riu) ifn'<u. 

By direct computation one verifies that the Green’s function satisfies the equation 

{H — A) sx{u, u') = 5{u — u') . 


(2.19) 


Thus taking s\{u,u') as the integral kernel of a corresponding operator s\ on TC^ = 
L^((0, tt), du), this operator is a bounded inverse of the operator {H — A). Thus A 
is in the resolvent set, and s\ is the resolvent. We conclude that the spectrum of H, 
defined as the complement of the resolvent set, is given as the set of all A for which the 
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Re V 



Im V 



Figure 1. Typical plot of the potential V . 


Wronskian (|2.18j) vanishes for non-trivial solutions satisfying the boundary 

conditions (I2.10h at u = 0 and u = tt, respectively. 


3. The Qualitative Behavior of the Spectrum 

We now explain qualitatively how the spectrum of the angular operator looks and 
how this qualitative behavior can be understood. This will also motivate and explain 
the statements in Theorem ll.il Before discussing the effect of the imaginary part, we 
consider the situation that and A are real, so that V is real-valued. Then the spec¬ 
trum can be understood most easily by considering the Sturm-Liouville equation (12.2h 
as a one-dimensional Schrodinger equation with Hamiltonian (12.3j) . As shown on the 
left of Figure [H the potential looks typically like a double-well potential. This po¬ 
tential is approximately symmetrie (because the quadratic terms in H are symmetric 
around 7r/2 according to (|2.5p . but the terms (12.61) are anti-symmetric). If instead 
of a double-well potential we had a single-well potential, the eigenvalue could 
be computed approximately for large n by the Bohr-Sommerfeld-Wilson quantization 
condition (see [13 §48] or [HI eq. (2.5.51)]) 

(b pdq = 27rn , 


where one integrates momentum along a closed classical path of the particle. Thus, 
denoting the zeros oiV = W — \ hy ui and Ur (with ui < Ur < Tr/2), we obtain 



W du = 2'Kn . 


(3.1) 


From this formula, the expected eigenvalue gaps can be computed by 


TT : 


j‘Ur(n+l) 

lui(n+l) 


rur{n) i-u, 

^Xn +1 -Wdu- V^n-Wdu^ / 

J U£{n) Jui 


2^\n-W 


du 


so that 


An-|-1 Af] 



1 

y /\ n-W 



-1 


(3.2) 


In particular, for the large eigenvalues we obtain Weyl’s asymptotics (see [211 Sec¬ 
tion 11.6]) 


An ~ , 


An-|-1 An — n 


(as n —>■ oo) . 


(3.3) 
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Re V Re V 




Figure 2. The poles of ReF in the case k = s (left) and k ^ s (right). 

Another spectral region of interest is if A lies near the minimum of the potential. 
Approximating VF by a quadratic potential and using that W ~ in this case we 
obtain the scaling 

XnC:in\Q\, A„+1 - An If^l (ifl<n<|n|). (3.4) 

These formulas describe the behavior of the eigenvalues if we had a single-well po¬ 
tential. The eigenvalues of the double-well potential can be understood by considering 
two Hamiltonians with a single-well potential and by weakly coupling them together 
via a potential barrier (see for example |2n( Section 3.3]). If our double-well potential 
was symmetric about vr/2, the two single-well Hamiltonians would have degenerate 
eigenvalues. Coupling them together slightly removes the degeneracy, leading to the 
well-known eigenfunctions with even and odd parity (similar as considered for example 
in m Sections 3.4 and 3.5]). In this way, we would end up with pairs of eigenvalues. 
These pairs would be separated by spectral gaps having the behavior (|3.2[) . Since in 
our situation, the double-well potential is not symmetric about vr/2, we do not know 
a-priori whether the eigenvalues of the two single-well Hamiltonians are degenerate or 
not. But we can conclude that the eigenvalues of the double-well Hamiltonian can 
appear at most in pairs, separated by gaps which again scale according to (j3.2[) . If A 
is chosen much larger than the potential barrier, the eigenfunctions no longer see a 
double-well potential. Therefore, Weyl’s asymptotics (13.3p should again hold for the 
large eigenvalues. 

These simple qualitative arguments already allow us to understand the statement 
of Theorem o in the special case of a real potential. Namely, the operator Qq is the 
spectral projection on all the small eigenvalues, for which the Born-Sommerfeld rule is 
not a good approximation. The operators Qi,Q 2 , ■ ■ ■ are spectral projection operators 
corresponding to one or two eigenvalues (depending on whether there is a spectral pair 
or not). 

Before moving on to the complex potential, we remark that in the case k = ±s, 
the potential at the pole goes to minus infinity (for a typical example see Figure [2]). 
However, it turns out that, using the known asymptotics of the wave functions near 
the pole, the above qualitative arguments still go through if we choose U£ r-^ |F!|-2 
close to the inflection point of the potential. 

We next discuss the situation for a complex potential. One potential method is 
to treat the imaginary part oi W as a slightly non-self adjoint perturbation (see m 
V.4.5] or as the application to the spheroidal wave operator in [8l Section 8]). For this 
method to be applicable, the imaginary part of the potential must be small compared 
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to the gaps, i.e. 


Imt^l < A„+i - An . 


(3.5) 


For any fixed hi, this condition will be satisfied for sufficiently large n in view of Weyl’s 
asymptotics (13.3p . But the inequality (|3.5I) cannot be satisfied uniformly in as the 
following argument shows: Using (|1.5p in (j2.4l) . one sees that 


sup I ImVFl > c |U| 

(0,7r) 


with a constant c which may be large. Therefore, the inequality (|3.5p is in general 
violated if we are in the asymptotic regime (13.4p . By choosing U large, one can arrange 
that this asymptotic regime includes arbitrarily many eigenvalues. We conclude that 
ImVF cannot in general be treated as a slightly non-selfadjoint perturbation. This 
means qualitatively that the imaginary part of W shifts the eigenvalues considerably 
on the scale of the gaps. The eigenvalues will typically move into the complex plane. 
Moreover, degeneracies and Jordan chains may form. 

In order to locate the spectrum in the complex plane, for a complex potential whose 
real part has a single well one can again use the Bohr-Sommerfeld condition (13.ip . 
which now makes a statement on both the real and imaginary parts of the integral on 
the left (|3.1I) . Treating the imaginary part of A — IT as a perturbation, we thus obtain 
to first order 


ru; 

J ue 

I 


y^Re (An — W) du = vrn 
Im(A„-IU) 


ui AyRe(An, — W) 


du = 0 . 


(3.6) 

(3.7) 


We will prove that these relations really make it possible to locate the spectrum in 
the complex plane. Applying these relations naively, we find for the Hamiltonian with 
the single-well potential that the real part of the eigenvalues behaves just as discussed 
for the real potential. The imaginary part of the potential, however, must be adjusted 
such that (j3.7l) holds. In particular, we again find that the spectral points form at 
most pairs, separated by spectral gaps which scale similar to (13.21) . If the two spectral 
points of the pair coincide, a Jordan chain of length at most two may form. In this 
way, one can understand all statements of Theorem 11.11 


4. Overview of the Proof of the Main Theorem 

Making the above qualitative arguments precise requires an intricate combination 
of different mathematical methods. In order to facilitate reading, we now give a short 
overview of the proof of Theorem 11.11 In Section [5l we collect general statements 
on Sturm-Liouville operators with a complex potential. We show that the spectrum 
is purely discrete, and that the Hilbert space can be decomposed into a direct sum 
of invariant subspaces. Moreover, idempotent operators mapping onto these invari¬ 
ant subspaces can be constructed using contour integral methods. In Section [6] we 
introduce a useful method for analyzing the oscillatory behavior of solutions of the 
Sturm-Liouville equation of the form ()2.17p . These estimates are essential for locating 
the spectrum and for making the Bohr-Sommerfeld condition (j3.6jl precise. 

Our proof involves a deformation argument where we continuously deform a real po¬ 
tential to our complex potential (Section [16]). Moreover, in our proof we will sometimes 
be able to treat the imaginary part of the potential as a perturbation (cf. Section [8|). 
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The starting point of these methods is to have detailed information on the spectrum 
and the spectral gaps for a real potential. These estimates are worked out in Section [71 

In Section [8] we employ the method of slightly non-selfadjoint perturbations to obtain 
the desired spectral representation provided that lies in bounded set (see Proposi¬ 
tion [8T]). Therefore, all the subsequent sections are devoted to the problem of getting 
estimates for large |n|, uniformly in the spectral parameter A. 

In Section [9| we derive an a-priori estimate for the imaginary parts of all eigenvalues. 
The method is to evaluate an expectation value (see (I9.ip l giving an equation which 
makes the Bohr-Sommerfeld condition (|3.7p precise. This a-priori estimate is needed 
in order to distinguish the different cases and regions in Section [TOl In Section [12] we 
shall return to the method and refine it considerably. 

For the remaining estimates we shall construct approximate solutions of the Sturm- 
Liouville equation by glueing together WKB, Airy and parabolic cylinder functions as 
well as asymptotic solutions near the poles at rt = 0 and u = tt. Moreover, we derive 
rigorous error bounds. In Section [TO] we give an overview of the different cases and 
regions and explain how to locate the spectrum. The detailed estimates are worked 
out in Section m 

Section gives refined integral estimates of the imaginary part of the potential 
(see Propositions 112. n and fT2l2p . These estimates make use of the specihc form of our 
potential and will be needed several times in the subsequent sections. 

In order to show that the Jordan chains have length at most two, our method 
is to show that if Aq is an eigenvalue, then there is an annular region around Aq 
which contains at most one other eigenvalue (see Figure [7]). In order to construct this 
annular region, we differentiate the equations with respect to A and use an implicit 
function argument. The A-derivatives are computed and estimated in Section [13] The 
construction of the annular regions is given in Section [T4l 

Section is devoted to estimates of the Green’s function. Here the main task is to 
estimate the Wronskian of the fundamental solutions in (j2.17p . 

In Section [16] we continuously deform the potential from a real potential to our 
complex potential. Combining all the results from the previous sections, we can track 
the eigenvalues and control the spectral gaps. We also derive uniform norm estimates 
for the operators Qn and show that their sum converges strongly to the identity. 

5. General Functional Analytic Results 

In order to get into the standard functional analytic framework, we consider the 
Sturm-Liouville operator (|2.3I) as an operator on the Hilbert space L^((0, tt)). As the 
dense domain of definition we choose those function in C'^((0, vr)) H L^((0, tt)) which 
satisfy the boundary conditions ()2.10l) . 

Lemma 5.1. The spectrum of the Hamiltonian ([23]) is discrete and has no limit 
points. 

Proof. For Sturm-Liouville equations with a continuous potential, this is proved in [U 
Chapter 12]. Since the potential (12.8|) has poles at u = 0 and u = n, we give the 
proof in detail. For any A E C, we choose non-trivial solutions 4>l and 4>n with the 
generic asymptotic behavior ()2.17p . These solutions can be chosen to depend locally 
holomorphically on A in the sense that every Aq E C has an open neighborhood U such 
that the functions (pL and (pn are holomorphic in A E 1/ (these holomorphic families 
can be constructed for example by taking the solutions of the Sturm-Liouville equation 
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for variable A but fixed boundary values at some u E (0,7r)). Then the fnnctions (/>£ 
and defined by (12.171) as well as their Wronskian in (I2.18P are also holomorphic 
in A E t/. 

Let ns show that the fnnction (j)^) does not vanish identically. If this were the 

case, by analytic continuation we would conclude that vanishes identically 

for all A E C. Thus for every A E C there would exist a non-trivial solution cj) satisfying 
the boundary conditions ()2.10p . On the other hand, the computation 

^|<^P = 2|<^f+ 2|<^|2 Re(ll^-A) 

shows that if A is large and negative, then the absolute square of 4> is convex away 
from small neighborhoods of the poles at 0, vr. But this convexity is incompatible with 
the asymptotics near the poles in (j2.10p . a contradiction. 

The result follows because holomorphic functions which do not vanish identically 
have isolated zeros. □ 


For a self-adjoint operator, one can construct the spectral projection operators by 
integrating the resolvent along a closed contour. In our non-selfadjoint setting, where 
the operator need not be diagonalizable, we cannot expect to obtain a spectral decom¬ 
position. But we can detect invariant subspaces: 

Lemma 5.2. Let T be a closed contour which lies entirely in the resolvent set and 
encloses points in the spectrum with winding number one. Then the contour integral 

Qr:=-^^'SAdA (5.1) 

defines a hounded linear operator whose image is the invariant subspace corresponding 
to the spectral points enclosed by T. 

The operator Qr is idempotent. Moreover, the product of two operators Qr and Qp' 
is given by 

QrQr' = Qf^ (5.2) 

where T is any contour which encloses precisely all the spectral points enclosed by T 
and r', all with winding number one. 


Proof. We first show that Qr is idempotent. Mnltiplying the identity 

{H- X)-{H - A') = A' - A 

for X X' from the left by sx and from the right by sy, one obtains the resolvent 
identity (see for example |18l Theorem VI.5]) 

^ X- X' ~ ' 


We let T' be a contour obtained by continuously deforming the contour T in the 
resolvent set snch that every point of T is enclosed by T' with winding nnmber one. 
Then, using the resolvent identity. 


Qr Qr 


Qr Qr' 



1 

A-A' 


(sa - Sy) 


1 



dX' A 
A-AV 


Sx dX -\- 





dX \ 
X-X'J 


sy dX' . (5.3) 
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Carrying out the inner contour integrals with residues, the integral in the first sum¬ 
mand gives —27ri, whereas the integral in the second summand vanishes. We conclude 
that Qr Qr = Qr- 

In order to prove the more general formula (15.21) . it is convenient to deform the 
contours and to decompose each contour integral into a finite sum of integrals where 
each contour encloses only one spectral point. Then, in view of the idempotence of 
the Qx, it remains to prove (|5.2I) in the case that T and T' enclose different points of 
the spectrum. By continuously deforming T' without crossing spectral points we can 
again arrange that the contours T and T' do not intersect. As in (|5.3p . we obtain 

Qr Qn = - jT ^ ,, dA + ^ rfV . (5,4) 

Since the contours enclose different points in the spectrum, no point of T is enclosed 
by r' and vice versa. Hence the inner integrals in (j5.4p vanish, proving that Qr Qr' = 0. 

It remains to show that the image of Qr consists of the invariant subspaces corre¬ 
sponding to all the spectral points enclosed in T. Let H C C be the open set enclosed 
by r. Since the spectral points are isolated, we may decompose H into a finite number 
of subsets such that the boundary of each subset is a closed contour enclosing only one 
spectral point. Thus it suffices to consider the situation that T encloses exactly one 
spectral point Aq. Since the Wronskian in (12.191) is holomorphic in A, the resolvent sx 
at Aq has a pole of finite order n. Iterating the identity 

HQr = (^(Asa - 1) dX = i£ XsxdX , 

Jy zttz Jy 

we obtain 

(H - Ao)”Qr = SA dA = 0 , 

where in the last step we used that the integrand is holomorphic. Hence every vector in 
the image of Qr is contained in the invariant subspace corresponding to Aq. Conversely, 
let '!/> be a vector in this invariant subspace. Then there is G N such that (iL — 
Xo)^'tp = 0, and thus 

0={{H-X) + {X- Xo)f^P = E (I) {H - X)H . 

fc=0 ^ ^ 

Multiplying by sx and solving for saV’j we obtain 

N ^ X 

= - E E ^ " Ao)-'= {H - Xf-^ V’ • 

k=l A A 

Taking the contour integral, a computation with residues yields 


QrV' 



{H - 
(A - Ao)" 


dX 



(A - Ao) 


dX 



N X X 

fc=o A A 






Hence ^|J really lies in the image of Qx- This concludes the proof. 


□ 
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This lemma also shows that the dimension of the invariant snbspace is at most the 
order of the pole of the resolvent. 

The next lemma bounds the resolvent away from the real axis. 

Lemma 5.3. If 

inf I Iml/| > 0 , (5.5) 

(0,7r) 

then the Wronskian in (|2.18l) has no zeros. Moreover, the resolvent is bounded by 

IIsaII < 1 inf I Im 1^1 
V(0'^) 

(where || • || denotes the sup-norm on Hk). 

Proof. Since ImV is continuous, the condition (15.511 implies that lm.V is either always 
positive or always negative. We only give the proof in the first case because the second 
case is similar. If the Wronskian in (|2.18l) is zero, there is a non-trivial solution (f of 
the Sturm-Liouville equation (12.711 with the asymptotics (I2.10p . In the case kPs 
differentiating the asymptotics (|2.11ll . one sees that 

(f>{u) ~ (l -|- 0(n)) and 4>'{u) ~ (l -)- 0(u)) . 

As a consequence, we do not get boundary terms when integrating by parts as follows, 

0 = {(j)\ {-dl + V)(f)L2 = {du(j) I <9n</')L2 + {4>\v(I)) l2 . (5.6) 

Taking the imaginary part, we conclude that 

PTT 

0 = / Imy|(^|2, (5.7) 

Jo 

in contradiction to (j5.5ll . In the case k = s, the situation is a bit more subtle because 
differentiating the asymptotics (I2.12|l . one sees that 

(j)(u) = Cl^/u + 0{u) and (f)'(u) = —^-|-0(u'^). 

2^/u 

This implies that integrating by parts in (15.6p we get real-valued boundary terms, so 
that (j5.7p again holds. The remaining case k = —s is treated similarly by differentiat¬ 
ing (EISI). 

Next, setting i/’ = s\(p and again integrating by parts, we obtain 
M\ U\\ > = \{'i(\i-dl + V)'il;)L 2 \ > |lm(V’|I^'!/>)i, 2 | > ||V^f inf Iml/, 

(0,7r) 

implying that 

\\(t>\\ > ||sa</>|| inf Iml/ . 

(0,7r) 

Since this inequality holds for all cf G (K^., the result follows. □ 

We are now in the position to state a general completeness result. The method is 
based on an idea in [H proof of Theorem 2.12] and was used previously in [9]. First, 
we write the potential (12. 8p and (|2.9p in the form 



F = IT- A 
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with W independent of A. For given i? > 0, we consider the two contours Fi and r 2 
in the complex A-plane defined by 

Fi =5SR(0)n|lmA < - inf \lmW\ - ^/R} 

I (0,7r) J 

r 2 = 55K(0)n|lmA> inf \ lmW\+VR}. 

I (0,7r) J 

and set T{R) = Fi U r 2 (see Figure [3|). 

Theorem 5.4. For any (j) G C'“((0, tt)), 

d>(u) = --—: lim 
27r'l R—>-oo 


/ ( 

Jr{R) 


isx(f))iu) dX . 


(5.8) 


Proof. Since the length of the contour U S 2 ■= dBii{0) \ T{R) only grows like x/R, 


r dX f 

JdBniO) ^ Jr(R) ^ R Js 


As a consequence, 


T(R) 


—: lim 
Sttz R —^00 


/ 

JV{R) 


R 

T 


S'iUS '2 


= 1 . 


|dA| 0 


(5.9) 


Since our contours lie in the resolvent set, we know that for every A G T{{R), 

f = sx{-dl + W-X)cj). 

Dividing by A and integrating over T{R), we can apply (|5.9I) to obtain 
^ iiSo X («(- al + tr - A) 0 ) (») 


-: lim 

27r2 R —^cxD 


f 

Jr{R) 


{s\4>)iu) - j{sx{-dl + W)(p){u)\ dX . 


But the second term in the curly brackets vanishes in the limit, because by Lemma [5.3l 

[ {sx{-dl + W)<p){u)^ < f 
Jt(r) ^ -'r(R) VR 


C |dA| ^ 2ttC 


|A| 


xfR 
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Figure 4. The osculating circle to the curve C{u). 


Thus (|5.8I) holds. □ 

This theorem shows that the operators Qr defined by ()5.ip converge to the identity if T 
tends to a contour which encloses the whole spectrum. The advantage of this method is 
that it does not require a functional analytic framework, but only uses properties of the 
Green’s fnnction s\{u,u'). The drawback is that one obtains strong convergence only 
on a the dense snbspace of test functions. In order to prove strong convergence on the 
whole Hilbert space, we will rely on the theory of slightly non-selfadjoint pertnrbations 
(see Section [5] and Section flfi.dp . 


6. An Osculating Circle to the C-Curve 


In order to locate the spectrum, we need to find the zeros of the Wronskian in (|2.18p . 
The main difficulty is to understand the behavior of the integrals in ()2.17p . To this 
goal, we now develop a method referred to as the “osculating circle method.” For ease 
in notation, we only consider the solution and omit the snbscript L. We denote 
the integral in (j2.17p by 

Then 

= ( 6 . 2 ) 
making it possible to relate the behavior of (jP to properties of the function Ci^)- 
In order to clarify the evolution of the function ('(u) in the complex plane, it is 
usefnl to consider the osculating circle to the curve at a point ({u) (see FignreHp. 
The curvature K of the curve C{u) and the radius R of the osculating circle are given 
by (see for example [U Theorem 5.1.6]) 


Im(C"C') 

IC'P 



The center p of the osculating circle is 

P = C 


K 1C' 


Moreover, we introduce the angle i? as the argument of C>, 




(6.3) 


(6.4) 

(6.5) 


Then 


2 
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SO that ()6.4I1 becomes 

(: = P + ^e-^^^. ( 6 . 6 ) 

In order to simplify the computations, we always choose the phase and normalization 
of (j) such that 

Imy(uo) = 1 (6.7) 

for some uq which will be specihed later. Moreover, we set 

2/0 = y{uo) ■ ( 6 . 8 ) 

Using the dehnition (16.ip of C as well as the differential equation (12.71) . we obtain 
useful formulas for i/, K, p and their derivatives. 

Lemma 6.1. 

i}' = lmy . (6.9) 

Furthermore, 


K{u) 

= 2\(j)\^ Imy , 

24,^ Im y 

(6.10) 

K'{u) 

= 2|(/.|2 ImU, 

p'(u) = —^ ImU . 

2 (fp Im^ y 

(6.11) 


R'{u) = 

ImU 

(6.12) 


2 (/>p Imy Imy 


Before coming to the proof, we point out that for a real potential, this lemma shows 
that K, p and R are constant. Thus for a real potential, the curve (^{u) lies on a fixed 
circle with radius R centered at p. The position of C on the circle is described by the 
angle t), and its evolution is described completely by (16.91) . Moreover, we point out 
that Ip'I = which can be understood from the fact that the point C stays on the 
circle as the osculating circles move. 


Proof of Lemma 16. IL A direct computation yields 








Id If” 

-Im log C' = -IiR — = Im 

2 du 2 C' 


Im(C"C') = -2 


Imy 



Imy 


This gives (j6.9l) as well as the formula for K in (|6.10p . Using this formula in (16.4p . we 
obtain the second equation in (|6.10l) . 

Next, we write out the real and imaginary parts of the Riccati equation (|2.14p . 

Re y' = Re U — Re^ y + Im^ y 
Im y' = Im U — 2 Re y Im y . 


Using ()6.14p . we obtain 

K' = 2|(/)p Im(y') + 4Rey \4>f Imy = 2 ImU , 


(6.13) 

(6.14) 
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giving the formula for K' in (|6.1ip . Moreover, 

P' = C' + 


* T / /N 1 i 

lm{y ) — + 


2 Im^ y 


(jp' 2 Im y (f)^ 


f Im y — 2 Re y Im y') — -|- 


— —T + 


— + 


2 Im^ y 
i Im V 


cjr 2 Im y (f)' 


To y 


cjP' 2<jP‘lTC?y (/i^Imy 
Finally, we differentiate (|6.3p . 

R' = - 

and using (jB.lOp and pB.lip gives (|6.12l) . 


* * 
Rey + 


ilmR 


(p'^lmy^ 2(lP‘lvc?y 


KK' 

Ts" ’ 
\K\2 


□ 


7. Estimates for a Real Potential 

In Section [TU] we shall consider a homotopy of the potential which joins the poten¬ 
tial V with its real part. In preparation for this analysis, we now derive eigenvalue 
estimates for a Sturm-Liouville equation with a real potential. More precisely, we 
replace the potential in the Sturm-Liouville equation (12.7p by its real part, 

where V is again given by (12.8p and (12.9p . We can assume that A is real, so that the 
equation can be written in the Schrodinger form (12.21) with the Hamiltonian 

(R 

H = -— + ReW . (7.2) 

du^ 

This Hamiltonian has a unique self-adjoint extension, as the following consideration 
shows: In the case k ^ s, the asymptotics in (|2.1ip shows that one of the fundamental 
solutions is square-integrable near u = 0, whereas the other fundamental solution is 
not. Using Weyl’s notion, the Sturm-Liouville operator is in the limit point case at u = 
0 (see [H Sections 9.2, 9.3]. In the case A: = s, on the other hand, according to (|2.12|) 
both fundamental solutions are square integrable. This is the so-called limiting circle 
case (see [H Sections 9.4]). In all of these cases, our boundary conditions (|2.10p give 
rise to a unique self-adjoint extension (for details see [Sl Sections 9.2, 9.3, 9.4] or |5l 
Chapter XHL2]). 

For ease in notation, we denote the selfadjoint extension of (|7.2[) again by H, and 
its domain of definition D{H). For the analysis of the spectrum, it is again useful to 
consider the Riccati equation corresponding to dni), which we write as 

y' = ReU-y2, (7.3) 

where we again set y := (j)'/4>. We consider complex-valued solutions of this equation. 
A direct computation (see also [S] eq. (3.8)]) shows the product \4>\‘^ Imy is a constant, 

w := \cj)\‘^ Imy = const. (7.4) 

This implies in particular that the function y cannot cross the real axis. 
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7.1. A Node Theorem. The classical node theorem (see for example [231 Theo¬ 
rem 14.10]) states that the eigenfunction of a Sturm-Liouville operator has ex¬ 
actly (n — 1) zeros. We now state and prove this node theorem in our setting. In the 
subsequent Sections 17.21 and 17.31 we will apply the node theorem to obtain eigenvalue 
estimates and the Weyl asymptotics. There are two reasons why we decided to give 
the proof of the node theorem in detail. First, due to our singular boundary condi¬ 
tions, the proof given in most textbooks does not apply to our problem. Second, our 
proof works with osculating circles and complex solutions of the corresponding Riccati 
equation. It can be used as an introduction to the methods needed later in this paper. 

Proposition 7.1. The spectrum of the Hamiltonian (|7.2p is a discrete subset of M. 
which is bounded from below. Numbering the eigenvalues in increasing order, Aq < 
Ai < ..., the eigenfunction corresponding to An has exactly n zeros on the open in¬ 
terval (OjVr). Moreover, choosing A = An, any solution y of the Riccati equation ()7.3I) 
with Im y > 0 satisfies the relation 


Imy = (n -|- 1) vr 


(7.5) 


Proof We let cfi and (j)2 be two real-valued fundamental solutions of the ODE (17.11) . 
Since their Wronskian 

w := - 4>l<p2 

is a non-zero constant, the functions 4>i and 4>2 cannot have common zeros. Hence 
the complex solution (p := 4>i + i4>2 has no zeros. By choosing suitable fundamental 
solutions, we can arrange that the corresponding solution of the Riccati equation (|7.3p 
satisfies (I7.4p with rc = 1, so that 

\(p\^ Imy = 1 . (7.6) 

The relations (|6.2p and (|6.ip define a solution cjp which satisfies the boundary 
condition at n = 0. The boundary conditions at u = tt are satisfied if and only 
if C('^) = 0. Hence the condition for an eigenvalue can be stated as 

C(vr) = 0 . (7.7) 

In order to control the behavior of the function <p, we again use the osculating circle 
method of Section [6l For a real potential, the relations ()6.1ip show that the center 
and the radius of the osculating circle are fixed. Moreover, combining the hrst identity 
in ()6.10p with (17.6p . one sees that K = 2. Hence the formula (16.6p simplifies to 


C(u) = p + ^ 


(7.8) 


where id satisfies the differential equation (16.91) . As a consequence, the eigenvalue 
condition (HZl) can be written as 


f 


Im y E ttZ . 


(7.9) 


The above formulas are valid if we let y be any solution of the Riccati equation in 
the upper half plane and if we satisfy (|6.7h (and consequently also (|7.6p ) by letting 


(p{u) = 


1 


Imy(Mo) 


exp 


UQ 


(7.10) 
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We now consider in particular a family of solutions y parametrized by A E M such that 

lim d\y{u) = 0 for all A E M . (7-11) 

U /^TT 


Such a family exists in view of the asymptotics near u = vr as worked out in m 
Section 8] (namely, one chooses (j) with the asymptotics as in [121 Section 8] with 
coefficients adjusted such that the leading asymptotics is independent of A, implying 
that d\y vanishes to leading order). 

Differentiating (|2.14l) and (16.8p with respect to A and using that dxV = —1 gives 


yA = -1 - 2yyA , y\{uo) = dxyo . 

Solving this linear ODE by integration, we obtain 

nu 

4 >^y\ = ^Ayo - , 

J Uo 

so that 


y\(u) = 


(p {uo) 


dxyo - 


1 n 

W)Juo 


-r v-y Juo 

Integrating this differential equation with respect to A yields 

ci: -/'(. 

In the last integral we perform the transformations 

PIT PU P7T P7\ 

/ du dv---= dv 

J un J ua J uc\ J V 


du 


' uq 


pUq PU puq puq puq pv 

/ du dv ■■■ = — du dv ■■■ = — dv 

Jo J uq Jo j u Jo Jo 


to obtain 


(7.12) 


(7.13) 


(7.14) 


du ■ ■ ■ 


( 1 /'U \ /‘TT PUQ 

m L '^) ^ L * ■ 

Taking the limit uo tt and using ()7.1ip . the relation (j7.14p simplifies to 

P7T PUQ PUQ 

/ dxy= (p'^C = / (p(p^ . 

Jo Jo Jo 

The representation for (fxp^ derived in Lemma 17.31 below shows that the function <p(p^ 
has a non-negative imaginary part, and that its imaginary part is even strictly positive 
on a set of positive measure. Therefore, 



9a Ihi y > 0 for all A E M , 

showing that for our family of functions y, the integral (j7.9p is indeed strictly increasing 
in A. 

Combining this strict monotonicity of the integral (j7.9p with the continuous de¬ 
pendence on the parameter A, the intermediate value theorem gives rise to eigen¬ 
functions \n which are uniquely characterized by their number of zeros. Since the 
integral ()7.9p is strictly positive, converges to zero as A —)• —oo and tends to infinity 
as A —oo (using the WKB asymptotics), we conclude that there is a sequence of 
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eigenvalues Aq < Ai < ... and that the eigenfunction corresponding to A^ has pre¬ 
cisely n zeros in the open interval (0,7r). Moreover, we conclude that (|7.5p holds for 
the family of functions y satisfying (|7.11l) . 

In order to show that (|7.5I) holds for any smooth family of solutions y with Imy > 0, 
we use the following continuity argument: For any fixed A = A,i, we denote the solution 
satisfying (I7.1ip by yo, and let y be any other solution with Imy > 0. For any r E [0,1], 
we let (yr) be the family of solutions of (I7.3p with initial conditions 

yr(f) = r2/o(f) + (l-T)?/(f). 

Then the condition (17.9p is satisfied for any r, 

Im yr E ttZ . 

By continuity, this integral is independent of r. We conclude that m holds for 
any r E [0,1], and in particular for r = 1. □ 



Remark 7.2. We remark that for a real potential, the eigenvalue condition (I7.9p can 
also be understood without going through the osculating circle estimates, as we now 
explain. Since the real and imaginary parts of (p form a fundamental system, the 
solution can be represented as 

(ff = c Im(e“*“(/)) (7.15) 


for a suitable phase a and a complex prefactor c. The zeros of (pP are then determined 
by the phase of (p, 

(p^ = 0 arg (p ^ a + ttZ . 

In particular, for (p^ to satisfy the Dirichlet boundary conditions, it follows that 


Differentiating gives 


— arg 
au 


argf^lg E ttZ . 


— Im log (p = Im — = Im y , 
au (p 


(7.16) 


and applying the fundamental theorem of calculus again gives dZS]). Moreover, one 
sees again that the integral in (|7.9p gives tt times the number of zeros on (0, vr) plus 
one. Using the osculating circle method has the advantage that with (j6.2l) we have an 
explicit formula for (p^, making it unnecessary to think about how the angle a in (I7.15|) 
is to be chosen. 0 


We append the lemma which shows that the integrand lm[<p(p^) has a definite sign. 

Lemma 7.3. For every solution cp satisfying the normalization condition (16.7p (for 
any uq E (0, tt)), 


1 


lni{(p{u)(p^{u)) = - \cp{u)\ 


1 - 


(p{u) 


lim 




ri)P •i;\0 (p{vY 


Proof. Clearly, <p^ is a linear combination of the fundamental solutions cp and (/>, i.e. 

(/>° = a(p + f3 (p (7.17) 

for suitable coefficients a, /3 E C. In order to compute these coefficients, we compute 
the Wronskians of (p^ with both (p and cp. First, using the ansatz (|7.17p . we get 

w{(p,cp^) = /3w{cp,'^) , w(^,4P) = -aw{(pf(p) . 
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Next, using the representation (|6.2I) and (|6.1I) . we obtain 

(t>^) = (p'cpC - (k {(PCY = -(fY C' = -1 

0 (00' = -'w{4>A)C- I0PC' = -w{4>A)C- • 

Comparing these formulas, we can compute a and /?. We obtain that for any v G (0, tt), 

a = Civ) + -, /3 =-^ . 

Finally, the normalization condition (lEZI) implies that t(;(0,0') = 2i. Computing 
Im00° using the above relations, taking the limit u 0 and using that hm„^o Civ) = 
0, we obtain the result. □ 

7.2. Lower Bounds for Small Eigenvalues. In order to obtain eigenvalue esti¬ 
mates, we need to count the number of zeros of the function (fP. Our method is to 
decompose the domain (0, tt) into subintervals on which the potential Re V has a def¬ 
inite sign. On every interval where ReR is positive, the number of zeros is a-priori 
bounded: 

Lemma 7.4. If I is a closed interval with ReRj/ > 0, then fP has at most one zero 
on I. 

Proof. Assume conversely that there is more than one zero on I. We choose two 
neighboring zeros ui < U 2 . Then, possibly by flipping the sign of 0 we can arrange 
that 0 °|(ui,u 2 ) >0- As a consequence, the function 0° is convex on \ui,U 2 \ (for details 
and other estimates using this convexity property see [U Section 5]). This implies that 
fp is non-positive on (ui,U 2 ), a contradiction. □ 


Lemma 7.5. Let y he any solution of the Riccati equation dZJ]) with \my > 0. Then 
the number Z of zeros of fP on an open interval I C (0, tt) is bounded by 

-1-1- [ Imy < Z < 1-|- [ Imy . 

JI JI 

Proof. One method of proof is to consider the osculating circle for a real potential (|7.8p 
and to note that the change of the phase i? is given by the differential equation (16.9p . 
Finally, the representation (|6.2I) shows that the zeros of fP coincide with the zeros of C- 
An alternative method is to use the representation (j7.15l) with 0 according to (j7.10p . 
and to make use of the fact that arg 0 satisfies the differential equation (I7.16P . □ 

Combining the node theorem of Proposition 17.II with the last two lemmas, we obtain 
the following corollary. 


Corollary 7.6. Let R,..., /^ C (0, vr) be open intervals such that Re V is non-negative 
on the complement o/ R U • • • U R. On the R we choose any solutions of the Riccati 
equation (1731) with Imy^ > 0. Then for the eigenvalue Xjy, 


TT iN — 2k — 1) < 


k „ 
i=l •'L 


Imyi < TT {N + k) 


We now apply this corollary to the spheroidal wave operator. We restrict attention 
to lower bounds for the eigenvalues, but remark that upper bounds could be derived 
with similar methods. 
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Proposition 7.7. For every constant C 3 > 0 and any parameters k,s, there is N = 
N{cs, /c, s) E N such that for all O in the range (|1.5I) with |0| sufficiently large, the 
eigenvalue is bounded from below by 

Xn > C3 |^^| . (7.18) 


Proof. In order to prove (|7.18p we consider A < C 3 |II| for large |II|. Then, due to the 
summand sin^ u in (j 2 . 8 p . the real part of the potential is non-negative except in a 
neighborhood of u = 0 and u = tt. By symmetry, it suffices to analyze the behavior in 
a neighborhood of m = 0. Then the estimate 

1 inP 

ReP > - (C 3 + l)\n\ + ^ + \n\^ (D{u^) (7.19) 

shows that the potential is positive if 

u > ui := 2 {c3 + 1 ) \n\-^ . (7.20) 

We begin with the case k s. In this case, the estimate (|7.19l) is improved to 

Re V > — (c3 + i)|n| + + |n|^ ©(n^). 

In particular, we conclude that for large |II|, 

ReR > -2(c3 + I)|0| {if k ^ s) . (7.21) 


We choose a (possibly empty) interval (u_,u+) such that ReR is negative inside and 
non-negative outside this interval. In view of (17.2011 we may choose < ui. In 

order to count the zeros of on the interval (tt_,n+) we can assume that the mini¬ 
mum uq of RelA lies in the interval (u_,u+) (because otherwise the interval (u_,u+) 
is empty, and there is nothing to do). We consider the solution y of the Riccati 
equation (17.3p with initial conditions 

y{uo) = iVl ReR(uo)| • (7.22) 


We now apply the T-method (see [121 Theorem 3.2] or [ 6 l Lemma 4.1]), choosing a = 0. 
Then U = ReV, a = 0 and 

_ ReW 
“ 2 

(see [I 2 I eqns. (3.3)-(3.5)]). Using that ReU is monotone increasing on [uo,u+) and 
monotone decreasing on (u_,uo], we obtain 


T{u) 


Re V (uq) 
Re V{u) 


1 

2 


giving rise to an invariant disk estimate with center m{u) = il3{u) and radius R given 

by 

f{u) = ^ (iReU(uo)l^ + |ReU(u)| iReU(uo)r^) 

R(u) = ^ (iReU(uo)l^ -|ReU(u)| iReU(uo)r^) ■ 

In particular, one sees that 


Imy < v^|ReU(Mo)| on (u-,U-f-) . 
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As a consequence, 

_ _ 

Imy < a/I Rey(Mo)| \u+ - U-\ < y/\ Rey(Mo)| ui < (2(03 + 1)) 


f 


where in the last step we applied (|7.20l) and (I7.2ip . 

In the case A; = s, on the other hand, the function Rel/ tends to minus inhnity 
as tt \ 0. We choose tt_ = 0 and < ni with ReI7(u+) = 0. We first apply the 
invariant disk estimate near the pole as worked out in [121 Section 8.1]. This estimate 
applies up to some u < Choosing u such that ReR is monotone increasing on 

the interval {u,ui), on the remaining interval [u,ui] we can again use the T-method 
with a = 0. Again, this gives rise to the estimate 



Imy < C(c3) . 


Working out similar estimates near u = vr, we can apply Corollary 17.61 with k = 2. 
We conclude that we can choose N such that for all sufficiently large |n|, the chosen A 
is smaller than the eigenvalue. This concludes the proof. □ 


7.3. Weyl’s Asymptotics. In the next lemma we show that our boundary condi¬ 
tions (I2.10|) give rise to the usual Weyl asymptotics. 

Lemma 7.8. The spectrum of the Hamiltonian (|7.2p with boundary conditions (|2.10|) 
lies on the real axis and consist of points Aq < Ai < • • • . For large n, the eigenvalues 
and gaps have the asymptotics 

Xn = n'^ + 0{n) 

An-i-i — An = 2n -|- O(n^) . 


Proof. We consider the family of solutions y of the Riccati equation (17.3p with initial 
conditions 

y(|) = i\/A for A € M"*" . 

Asymptotically for large A, the potential ReR becomes nearly constant according 
to (12.811 . except at the poles at u = 0 and u = n. In the case k ^ s, one can control 
the behavior near the poles by using the T-method similar as explained after (I7.22p . 
We thus obtain asymptotically 


rn 

/ Im y = a/A vr -|- 0 

Jo 


1 

7x 


(7.23) 


In the case k = ±s, one can use the asymptotics of the fundamental solutions as 
worked out in m Sections 7 and 8] to again obtain ()7.23p . 

Combining (|7.23p with (I7.5p gives the result. □ 


8. Slightly Non-Selfadjoint Perturbations 

We now prove Theorem 11.11 under the additional assumption that 11 is restricted to 
a bounded set; 

Proposition 8.1. Let U G C he a bounded set. Then for any s and k in the 
range m, there is a positive integer N and a family of operators QniLl) on “Kk 
defined for all n G'NL) {0} and Ll G U which has the properties (i)-(v) in the statement 
of Theorem \l.l[ 
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Figure 5. Contour integrals for slightly non-selfadjoint perturbations. 


This proposition differs from Theorem 11.11 by the fact that here the parameter N may 
depend on the set U, whereas in Theorem ll.li the parameter N is to be chosen uniformly 
for all n in the unbounded strip (jl.5|] . This uniformity in is the main difficulty of 
the present paper; its proof will be the concern of the remaining Sections fTBlllhl 


Proof of Proposition \8.1{ We again consider the Hamiltonian (17.21) with a real poten¬ 
tial with boundary conditions (j2.10l) Choosing contours which enclose each of the 
eigenvalues Ao,Ai,... with winding number one, the contour integral m dehnes 
idempotent operators Qn, n G 0,1,.... Since H is formally self-adjoint, these oper¬ 
ators are symmetric, implying that the Qn are orthogonal projection operators onto 
mutually orthogonal subspaces. For any A G C \ {Aq, Ai, ...}, we dehne the resolvent 
of the self-adjoint problem by 

oo ^ 

= E ■ 

_n ''' 


Here the sum converges absolutely in L(fH). Moreover, the resolvent satishes the 
identities 


du^ 


-|- Re IT — Al ) sa = 1 


and IIQaII = sup , ^ 

neNU{0} 1^ 


Our method for treating the imaginary part of the potential is to use the theory 
of slightly self-adjoint perturbations (see [151 V.4.5]), similar as worked out in [8l 
Section 8] or |3l Chapter 12]. We hrst note that, since the poles in ()2.5I) are real¬ 
valued, the imaginary part of the potential is bounded. 


|ImlT(n)|<C' for all n G (0, tt) and H G H (8.1) 

(where the constant C clearly depends on U). Next, using Weyl’s asymptotics of 
Lemma 17.81 we can choose N so large that 


|A„+i — A„| > 4(7 for all n > . 

We choose contours T^ (for n > N) as circles centered at \n with radius 2(7. More¬ 
over, we choose Tq as a circle which encloses the eigenvalues Aq, ..., Aw-i, and whose 
distance to the spectrum is at least 2(7 (see Figure [5]). Then for any A on one of these 
contours, 


1 


( 8 . 2 ) 
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This makes it possible to define the resolvent for the Hamiltonian (12.311 with the com¬ 
plex potential, which we denote for clarity by a tilde, via a Neumann series, 

CXD 

SA := ^ ( - SA IniTT)^SA . (8.3) 

k=0 

We now integrate this resolvent along the contours T^, 

Qn ■= ~ 2^ s\ dX , n G {0, N, N -\- 1, ■ ■ ■} ■ (8-4) 

As explained in [151 V.4.5], these operators are idempotent and map onto the invariant 
subspaces corresponding to the spectral points enclosed by the contour. Moreover, it is 
shown in [151 V.4.5] that the spectral projections are complete. The bound (II.6p follows 
immediately by estimating the contour integral (18.4h and the Neumann series (18.3h 
using (|8.2p and (|8.ip . □ 


9. An A-Priori Estimate for Im V 

Assume that A G C is an eigenvalue. We let (JP be a corresponding eigenfunction. 
This function satisfies the Dirichlet boundary conditions at u = 0 and u = tt. There¬ 
fore, the corresponding functions cj)^ and as defined by (I2.17P are both multiples 
of (jP. Multiplying the differential equation for by (jP and integrating, we obtain 

»= (- 1 ? + • <"> 

where in (a) we integrated by parts and used the asymptotics for the decaying solution 
in (12.lip and (I2.12p . (I2.13p to conclude that the boundary terms vanish. We thus obtain 
the relation 

ImV|(/.°|2 = 0. (9.2) 

This identity immediately gives rise to the following a-priori estimate. 

Lemma 9.1. Suppose that X £ C is an eigenvalue. Then 

|ImA|, |ImV| < C|H| 

with a constant C which is independent of X and Q. 

Proof. Using the explicit form of the potential (|2.8p together with ()1.5I) . one sees that 

I Im V — ImA| < 2c|H|-(-const . (9.3) 

The integral equation (19.21) implies that the function ImV must change sign on the 
interval (0,7r). As a consequence, the absolute value of ImA is bounded by the right 
side of (19.3p . This gives the result. □ 
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10. Overview of the Estimates for a Complex Potential 

We now enter the general estimates. Recall that our equations involve the parame¬ 
ters k, s, Q and A. We always keep k and s fixed. The parameters Q and A, however, 
may vary in a certain parameter range to be specified later on, and we must make 
sure that our estimates are uniform in these parameters. In order to keep track of the 
dependence on Q and A, we adopt the convention that 

all constants are independent of and A , 

but they may depend on k and s. Moreover, in order to have a compact and clear 
notation, we always denote constants which may be increased during our construc¬ 
tions by capital letters Si, 621 • • •• However, constants with small letters ci, C 2 ,... are 
determined at the beginning and are fixed throughout. We use the symbol 

< • • • for < c • • • 

with a constant c which is independent of the capital constants 6 / (and may thus be 
fixed right away, without the need to increase it later on). 

When increasing the constants C;, we must keep track of the mutual dependences 
of these constants. To this end, we adopt the convention that the constant S; may 
depend on all previous constants Ci,... ,C/_i, but is independent of the subsequent 
constants Cj+i,.... In particular, we may choose the capital constants such that Ci <C 
62 • • • . This dependence of the constants implies that increasing C; may also 

make it necessary to increase the subsequent constants , C/+ 2 , • • •• For brevity, 
when we write “possibly after increasing C;” we implicitly mean that the subsequent 
constants C^+i, C;+ 2 ,... are also suitably increased. 


10.1. Different Cases and Regions. In view of Proposition 18.II it suffices to con¬ 
sider the case that |H| is large. Thus in what follows we always assume that 

|H| > 64 . 


Since the imaginary part of H is bounded by (II.5p . by increasing C 4 we can always 
arrange that 


ReQf > ^ |H|2 . 


Furthermore, Lemma 19.11 gives us an a-priori bound on the imaginary part of the 
eigenvalues. 


ImA<|L!|. 


( 10 . 1 ) 


Moreover, in view of Proposition 17.7( we know in the case of a real potential that by 
choosing N sufficiently large, it suffices to consider the case that A is real and A S> |H|. 
With this in mind, in our estimates we may restrict attention to the case 


ReA > CsIHI . 


( 10 . 2 ) 


This inequality will be justified a-posteriori by showing that if we deform the potential 
continuously starting from a real potential and ending with our complex potential V, 
then the inequality (II0.2p will be preserved for all spectral points Xn with n > N (for 
details see Section flb.Sp . 
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For large |n|, the real part of the potential looks qualitatively like a double-well 
potential (see the left of Figured]). More quantitatively, in the region [f,^] away 
from the poles at n = 0 and rt = vr, according to (I2.8p we have 

Rel/ = Re(02) sin2u-ReA + 0(0) (10.3) 

Re R'= 2 Re(n^) sinrt cos n-|-0(0) (10-4) 

Re V” = 2Re(02) cos(2m) + 0(0) . (10.5) 


In particular, one sees that Re V has a unique local maximum at a point near which 
we denote by Umax, 


ReV'(umax) = 0 and n^ax = ^ +0(1^1 • 

Moreover, the real part of the potential is concave near this maximum, 

|Op 


- 2|0|^ < ReR" < -- 


on 


TT 277 

L3’TJ 


( 10 . 6 ) 


As the intervals (0, Umax] and [umaxj tt) can be treated similarly, we mainly restrict 
attention to the interval (0,Umax]- The value of the real part of the potential at its 
local maximum distinguishes different cases: 


' WKB case if Re R(umax) < -Ci |0| 

< parabolic cylinder case if —Si |0| < ReR(umax) < Ci |0| (lOT) 

Airy case if ReR(umax) > Si |0| . 


Here Si is a new constant which later on we will choose sufficiently large. 

In each of the above cases, we estimate the solution by considering different regions, 
as we now explain. First, we distinguish the pole region as the interval (0, u^) with 


ue := 


Si 

VRc a 


( 10 . 8 ) 


To the right of the pole region, there is a (possibly empty) WKB region {ui,Ur). The 
definition of Ur depends on the different cases. In the WKB case, we simply set Ur = 
U ma x- In the parabolic cylinder case, the fact that the function ReR is concave (|10.6I) 
implies that there is a unique point Ur G (f ,Umax) with ReR(ur.) = —Si |n|. In the 
Airy case, we make use of the following result. 


Lemma 10.1. In the Airy case, there are unique points Ur,u^ in the interval 

u<u< min [umax,u) (10.9) 


with 

Re R {ur) = —n , Re R(u+) = n 

where u, u and v are defined by 


u = 


VRe A 

YfoT ’ 


u = A 


VRe A 


V = min (j(S 5 ReA)R ;^(Si ReR(Mmax))0 • 


( 10 . 10 ) 


( 10 . 11 ) 

( 10 . 12 ) 
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Proof. We first show that Rei^ is strictly increasing on the interval (u, rimax)- First, 
according to (|10.4I) , by increasing 64 we can arrange that the function Re V is strictly 
increasing on the interval (f, Moreover, the concavity of ReR implies that ReR 


is also monotone increasing on the interval [^, Un 
we have the estimate 
|2 


ReR' > 


> 


2 
1 
2 


u - 2 > 


w 






u 


16 c|f 2 p 

Re^A 




16 c 


On the remaining interval (u, ■ 


lOpM > > 0, (10.13) 


where in the last step we possibly increased 65 . We conclude that ReR is strictly 
increasing on the whole interval (n, Wmax)- 
Next, at u we have the estimate 


ReR(n) < ^ -ReA = 


Re A 4c|ll|^ 




^ ReA 
Re A 


— Re A 


Re A 4c ^ , 

<-hwn ReA — ReA< — 

- 4 e| - 2 ‘ 

where in the last step we possibly again increased 65 . Moreover, iiu<u„ 
the estimate 


we have 


1 


ReR(u) > 7 \n\‘^u^ - ^ - ReA > 4 ReA- 




4c|0|^ 

Re A 


— Re A 


4c 


> 4 Re A — Re A — Re A > Re A , 

where in the last step we possibly again increased C 5 . Next, it follows from the defini¬ 
tion of v in (jl0.12|) and (110.2^ that v < Re A/4. Hence 

ReR(u)<—and Re R( min (u, Umax ))> u . 

Now the existence of solutions Ur and u+ of (jlO.lOp follows from the intermediate value 
theorem. Uniqueness is an immediate consequence of the above strict monotonicity 
of Re R. □ 


To summarize, the point Ur is defined by 

Ur = Umax ln the WKB case 

< ReR(ur) = —Cl |H| in the parabolic cylinder case (10.14) 

Re R (ur) = —v in the Airy case . 

In the Airy case, the interval (u+,Umax) with u+ as in (jlO.lOp is another WKB region 
to the right of the zero of Re R. We thus obtain the following regions: 


pole region 

(0,U£) 

in all cases 

WKB region 

{U£, Ur) 

in all cases 

parabolic cylinder region 

(Ur , rZmax) 

in the parabolic cylinder case 

Airy region 

(Ur,U+) 

in the Airy case 

WKB region with Re R > 0 

(U-|-, Umax) 

in the Airy case . 


(10.15) 

The different cases are illustrated in Figure [6l 
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10.2. Locating the Eigenvalues. Our general strategy is to construct a special solu¬ 
tion yL of the Riccati eqnation (I2.14p on the interval (0, tfmax]) and a special solution yn 
on the interval [nmaxj'r). These solutions are defined by the initial conditions 

yL{uo) = Vo , Vr{uq) = yQ , (10.16) 

where Uq and Uq are chosen near the poles at tt = 0 respectively u = tt (for details 
see Section [11.21 and 111.31 below). We choose these special solutions in snch a way 
that our estimates become as simple as possible. This means in particnlar that these 
solutions have no singularities. Then we introdnce corresponding smooth solutions of 
the Sturm-Liouville eqnation (|2.7I) by integration (cf. (I2.15p ). 


(t)L{u) := exp 




(10.17) 


both normalized according to (IS21). These solutions will not satisfy the Dirichlet 
bonndary conditions (|2.10p . By introdncing the fnnctions (j)^ and again by (|2.17p , 
we obtain solntions which do satisfy the Dirichlet bonndary conditions. In order to 
locate the eigenvalues, we must analyze the eigenvalne condition (|2.18l) . It is most 
convenient to evalnate the Wronskian at ttmax) 




Similar to (|6.ip and (|6.2p we set 



4>l = 4>lCl, 


4>r = 4>rCr- 


(10.18) 


(10.19) 

( 10 . 20 ) 


Differentiating these relations, we obtain 
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and thus 


w{(I)lAr) = (t>L <pR 


D^/ 


i4>LJ 


D\/ 










= (t>L Cl (t>R Cr 2/L + 


VR 


1 


C>ICr 


VR 


1 


C^ICl <^^Cr. 

Therefore, the eigenvalue condition (jlO.lSp can be written alternatively as 

1 1 


{VL - Vr) + 


= 0 . 


( 10 . 21 ) 


( 10 . 22 ) 


C^ICl c>ICr 

Indeed, in this form the eigenvalue condition is most suited for our analysis. Our 
main task is to analyze the behavior of the functions yi and yR as well as the derived 
functions (l)L,4’R and ClXr (obtained by (I10.17p and (110.1911 . (I10.20|l i. 

11. Estimates in Different Regions 

11.1. Estimates in the WKB Region (ui,Ur). The name “WKB region” suggests 
that in these regions the WKB solutions should be a good approximation. This really 
is the case, in the following sense: 

Proposition 11.1. For any (5 > 0 and for sufficiently large Ci, the WKB conditions 

lE'l \V''\ \V"'\ 

[y, ^<<5 ( 11 . 1 ) 

|E|i |E|t 

hold in the WKB regions {u£,Ur) and (n+,nmax) (see (110.1511 ). uniformly in Q, and A. 

For the significance of the inequalities (lll.ljl we also refer to |12( eq. (4.1)] and the 
estimates in m Section 4]. 

The proof of this proposition is split up into several lemmas. The proof will be 
completed at the end of this section. 

Lemma 11.2. Possibly by increasing Ci, we can arrange that in the region u > ug, 
the potential and its derivatives are bounded by 


jwi < \n\^u^ + 


Re A 
2 

(Re A)2 




(ReA)i 

ef 


ef 


|E" 




< 


Proof. Using the explicit form of the potentials in (|2.5I1 . (12.611 and (|2.8p . we obtain 

|1E| < + 4 ^ + i ReA < + ^ 


el 


ju'l < 2\n\'^u + 




u. 


ef 




Up 


ef 






ef 


where in (a) we possibly increased e 


1- 


□ 
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Lemma 11.3. The statement of Provosition [7l.ll holds in the region 

U£ < u < u, 

where u is again defined by (jlO.lljl . 

Proof. Using (j2.8p . we obtain 


Hence 


|U(u)| > I ReA| - |tU| > ^ - |J1|V > ^^ > Re A 


( 11 . 2 ) 


(11.3) 


IV^'I < < |12| Vr^+ 


|U' 


^ < 
o rv_/ 


ef 

\n\ 11 1 


ef 


|y|t ReA ' Cf ~ 65 Cf 

l^"l < 


|Q |2 11 1 

+ 7^ < 7VT + 


|u |2 ~ (ReA )2 ef ~ e| ef 


\V" 


< 




11 1 

_I-<-1- 

5 ^ o5 ~ 5 1 ' 1 


|V'|i~(ReA)§ e 5 ~e|e| ^5 


giving the result. 


□ 


It remains to consider the complement of the region ()11.2p . This complement is 
empty unless 

\/Re A vr 

2\n\ - 2 ■ 

Therefore, in what follows we can assume that 

ReA<|b!| 2 . 

Moreover, 


u > 


VRe A ^ Y^Re A| ^ 


> 


> 


(11.4) 

(11.5) 


2|f2| ' 2|I2| 2v1^ ~ ■ 

Lemma 11.4. Under the assumptions of Proposition [71. il 

\V"\ \V"'\ 

Proof. From piO.lOp . (110.12h . (llO.lh and (I10.2P we know that 

|U| >i/>min(e5|f2|,ei|I2|) = Ci |f2| . 

Combining this inequality with (111.41) and (lll.5p . the derivatives of the potential can 
be estimated by 


rill\ ^ io |2 , f 


|U'"|<|f 2|2 + 


U'^ 


| 0|2 + | 0|2 <| 0 | 2 . 


As a consequence. 


W" 


< 




|U" 


|U|2~e2|f2|2 02’ 


< 






|Q |2 


5 ’ 
2 
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completing the proof. 

It remains to estimate the term involving the first derivatives in (|ll.ip . 
Lemma 11.5. The statement of Provo sition Ul.T holds in the region 

u>u, 

where u is again defined by (| 10 . 1 ip . 

Proof. We can bound the potential from below by 


□ 


ReV > ^-^u‘^-ReX > 

- 4 ~ 8 


■ u 


Hence 


W| ^ Ifipu ^ (Re A)2 


iHii ~ 


+ 


< 


|fl| 1 1 1 


ef 

giving the result. □ 

It remains to estimate the term involving the first derivatives in (|Il.ip in the region 

u < u <u . ( 11 - 6 ) 

We begin with a lemma in the Airy case. 

Lemma 11.6. The statement of Provosition Ul.T holds in the Airy case if 
(ei|II|2 ReA)^ <2(e?|fi|2ReR(u^ax))^. 

Proof. In view of pi0.12p . the assumptions imply that 


1 


z2=-(ei|0|2ReA)3. 


(11.7) 

( 11 . 8 ) 


Moreover, according to (jlO.lOp . we know that |17| > n. Hence 

1 


|V^'l<l^^l"n + 4<|0|2n(l + 

\n\^ 




1 + 


~ " ' Re^A 


ITMl i 

< 

rs_/ 


ao.iiii 


\n\^u < |o|Vr^ 


|W| ^ |ii| Vr^ ™ 1 


This concludes the proof. 


□ 


Thus in the Airy case, in what follows we may assume that (lll.7|) is violated. Since 
in the WKB region under consideration, we stay away from the zeros of Re V in the 
sense that | RelA| > n with n as in (jl0.12p . it follows that 


ReR| > Rey(u„,.)) 


in the Airy case 


(11.9) 


We now return to the analysis of the region (|11.6p . without specifying whether we 
are in the WKB, the Airy or the parabolic cylinder case. As a consequence of (lll.bp . 

\v'\<\n\^u + ^<\n\V^<\n\^, 
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where in the last step we again applied (111.411 . From this inequality, we obtain the 
desired estimate provided that one of the following two inequalities holds: 

|0|2 


|y(M)| > 02 |0|3 or Rey(Umax)> 


Namely, the first inequality implies that 

IR'I . 


< 


|R|2 e| 


1 

3" 

el 


( 11 . 10 ) 


( 11 . 11 ) 


On the other hand, if the second inequality in (| 11 . 10 ll holds, we are in the Airy case 
(possibly after increasing 64 ), so that (|11.9I1 yields the estimate 

|F| > (e2|0|2Rey(n„,,))5 >ef |0|t. 

1 . 

This implies that the first inequality in (jll.lOh again holds (for 62 = 6 ®), making it 
possible to again use the estimate (lll.lljl . 

It remains to consider the case that both inequalities in (111.1011 are violated, i.e. 


|F(u)| < 62 |n |3 and ReR(umax)< 


Cl 


In this case, | Re R(tt)—Re R(umax)| ^ |Op/Ci, implying that u ss Umax (more precisely, 
by increasing 64 we can make |u — UmaxI arbitrarily small, uniformly in Q and A). 
Since ReR is concave near Umax (I10.0|l - we may integrate this inequality to obtain 

inP 

\U - UmaxI < \V'{u)\ < |Op |u - Umax] 

(U - Umax)^ < Re R(Umax) “ Rc R(u) < (u - Umax)^ • 

Hence 

\V'{u)\ < |0| v'ReR(Umax)-ReR(u) < |0| (^1 Re R(«max)| + \/|R(i/)|) 

In the case | ReR(umax)| < |R(^^)I) d follows that 


|R'(u)| ^ 2|F!| 


< 


|R(u)|i \V{u)\ 


< 


1^1 

ei|F!| 


1 

^ ’ 


giving the result. In the remaining case | ReR(umax)| > |R(ii)|) we know from (I10.14h 
and (llO.Th that we are in the Airy case. Hence, again using (lll.9h . we obtain 


|R^(n)| ^ 2|0| v^ReR(ua 


< 

rs_/ 


22 |Q|y/|ReR(Umax)| 


2 ! 


|R(u)|t \Viu)\i |n|2|ReR(Umax)| Cl ■ 

This concludes the proof of Proposition 111.11 

11.2. Estimates in the Pole Region in the Case fc = s. In this section, we 
analyze the pole region (0,U£) in the case k = s. We consider the solution cj)]^ as 
dehned by (110.1611 and (|10.17p . For ease in notation, we omit all subscripts L. The 
parameter uq in (|10.16ll is chosen as 

uo := ^ \n\-^ . 


( 11 . 12 ) 
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Possibly by increasing 65 , we can arrange that uq > (cf. (llO.Sp and (|10.2I) '). so 
that no lies in the WKB region. We choose the initial values at uq in (110.161) equal to 
the value (/'wkb/'^wkb for the WKB approximation, 

yo = lm^/V{uo) > 0. (11.13) 

We expand the potential near n = 0, 

V{u) = - n + + 0{\n\u^) + 0 (| 0 | 2 n^) . (11.14) 

Lemma 11.7. For any (5 > 0, we can arrange by increasing C 5 that 



uniformly in H and A. 

Proof. Using the asymptotics as worked out in m Section 7.1], on the interval (0, no) 
the solution (f has the form 

4>{u) « -c^/n (Ko{y/JIu) + {argy/JI- log(2) + 7 + z) Io{y/J2u)J , (11.15) 

where c is the constant (see [HI eqn. (7.3)]) 

c = - 

arg^/^-log( 2 ) + 7 + Z 

(and n is related to A and Q by (I2.9p ]). As specified in [HI Section 8.1], the error 
in ()11.15p becomes arbitrarily small for large |U|. Note that c is bounded uniformly 
in fi. For small n, the function cj) has the asymptotics (see m Section 7.1]) 

= —c[y/u log\^u\ +iy/u + . (11.16) 

In particular, using that 

[ ^ N = arctan (log(V^n)) + | , 

Jo n(l + logU^n) 2 

one sees that 1/|</>P is integrable. Hence 

fV ^ 

lim / -—— = 0 , uniformly in u . 

y\oJo W 

In view of (|10.8p . by increasing 65 we can make y/Jiui as small as we like. This gives 
the result. □ 

11.3. Estimates in the Pole Region in the Case k s. We now analyze the pole 
region (0, ut) in the case k s. We again consider the solution (pL as defined by (|10.16[) 
and (I10.17h and omit all subscripts L. The parameter uq in (jl0.16p is chosen as the 
minimum of Re V, 

ReU'(uo) = 0. (11-17) 

Introducing the abbreviation 

A = (fe - s)2 - i , 

the potential near u = 0 has the expansion 

u(u) = — A + + o(|n|tt^) + o(|np . 


(11.18) 
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Computing the zero of the derivative, we obtain 

uo = A^|f^r5+0(|Ori), (11.19) 

so that for large we have the estimates 

^ A3 < 1(0 < 2Ai |Q|"5 . (11.20) 

Possibly by increasing the constant 65 in ()10.2p . we can again arrange that uq > u^, so 
that uq lies in the WKB region. We again choose the initial values at uq in agreement 
with the WKB approximation ()11.13p . 


Lemma 11.8. For any 5 > 0, we can arrange by increasing C 5 that 



uniformly in 11 and A. 


Proof. Using the asymptotics as worked out in m Section 7.2], the solution cj) has the 
form 

4>{u) « c^/uK\k-s\ ( - y/pu) , ( 11 - 21 ) 

where c is the constant 


c = \j - 

vr 


As specified in m Section 8.2], the error in (|11.21l) becomes arbitrarily small for 
large jUj. For small u, the function cj) has the asymptotics (see [T2l Section 7.2]) 


= c ■ 


(n- 1 ) ! f 
y/2 fii 


(1 + 0{u)) 


— c 


y /2 (2 


jfc — sj! /i4 


\ + \k—s\ 


( 11 . 22 ) 


(l + 0 (u)). 


In particular, one sees that ]())] has a pole at u = 0 and is thus bounded from below 
near u = 0. Hence 

r 1 

lim / ——K = 0 , uniformly in a . 

^\oJo W 

In view of (|10.8p . by increasing C 5 we can make y/Jiui as small as we like. This gives 
the result. □ 


11.4. Estimates in the Parabolic Cylinder Region. In the next proposition we 
estimate the Riccati solution in the parabolic cylinder region. 

Proposition 11.9. Assume that in the parabolic cylinder region [ur,Umax]? one of the 
following two conditions hold: 

(a) The potential has a positive imaginary part, Im> 0. 

(b) The imaginary part of the potential has a zero on [u^ ,u^]. 
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Moreover, assume that the Riecati solution begins in the upper half plane, lmy{ur) > 0. 
Then there is a constant C 2 (depending on Qi) such that for large |0|, the solution on 
the interval [tirj'Wmax] can he estimated in terms of y{ur) by 

\y{u)\ < 62 \yiur)\ 

lmy{ur) 

lmy{u) > --- 

*-2 

< \(t){u)\ < 62 \(l)iur)\ . 

Proof. We set v = Ci|r2|. Using that the function ReU is concave near Umaxi we obtain 

(Wmax - Urf < (11.23) 

and thus 

l' (ttmax - Urf < Sj ■ 

Our strategy is to estimate y using the T-method as introduced in [l2l Section 3.2] 
choosing 

a = and ^ = 0. (11.24) 

Hence 

V = = 21^ and U = ReU — . (11.25) 

In case (a), our method is to apply [T2l Theorem 3.3] for g = 0. The terms Ei,... 
are estimated as follows, 


E 2 — E^ — 0 

li^il <^|ReU-ReU| + ^<^^+^ 

/f V V 


l^sl 


^ |ImU| W \n\ 


< 

rsj 


where in (*) we used (|1U.1|) as well as the fact that 

im 

|lm(02)| <|0 ||lmO| < |0 |. 
As a consequence, we can apply Lemma lll.lOl to obtain 

CMmax ^ s 

Jur ^ 

1 
V 


Wmax 1^1 


< \/z^(^max — Ur) + — 

— \/z^ (Umax Ur) T ^ Re U(Uniax) R^ U(^^)^ T (^max Ur) 


| 0 | 


<61 + 2 + 61^ <61 + 3. 

V 


This concludes the proof in case (a). 

In case (b), the imaginary part of V could be negative. Therefore, in order to 
apply |12l Theorem 3.3] we need to choose the function g positive in accordance with 
the inequality 

5>r-l. (11.26) 
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We choose a and P as in (I11.24P and g = |ri| 2 . Then the error terms Ei, E 2 and £'3 
estimated just as above. Estimating Im V with the help of the mean value theorem by 

I Iml/| < (u+ — u~) sup \V'\ < |n| (uf — n^) 

(where in the last step we used the explicit form of the potential (I2.8p l. the error 
term £ 4 , is estimated by 


/ l^ 4 |=W 

t/ Zlr J Ur 


\lmV\ 


^ 9 


r'^mayi 
J Ur 


|0| {u^ — 


< 


9\^\ ( R _ L^2 9_^ ^ gV^l ^ 

^ ^ v1^ |b!|r 


This can be made arbitrarily small by increasing |n|, implying that the inequal¬ 
ity (I11.26P holds. This concludes the proof. □ 

11.5. Estimates in the Airy Region. We proceed with estimates in the Airy region. 
We first recall that it remains to consider the interval (jlO.bp . Eor this interval to be 
non-empty, we can again assume that (|11.4p holds, 

ReA<|n|^ and thus |ReE|<|np. (11.27) 

We begin with a preparatory lemma. 

Lemma 11.10. In the Airy region, the function ReE is strictly monotone. Moreover, 

(u+— ^ 64 . (11.28) 

Proof. The strict monotonicity was already shown in the proof of Lemma 110.11 In 
preparation for the estimate (|11.28l) , we recall that the region (ur,u+) is contained in 
the interval (n, u) (see Lemma flO.ip . and thus 

\/Re A , , 

Ur,U+7^ 1^1 ■ (11.29) 

We consider the regions (u, ^) and [^,Umax) separately. In the region {u, ^), we 
know from (llO.lSp that ReE' > Moreover, the method in (|10.13p also gives the 

reverse inequality, 

R.V<|qG + A<(i + ^)|n|A 


1 + 


Re^A 


< |Opti 


where in the last step we used (jll.29p and (llO.lh . We conclude that 

^ ~ VReA . 

As a consequence, the mean value theorems 

(u+ — Ur) inf ReE' <2v < {u+ — Ur) snp ReE' 
K-«+) K,n+) 

give rise to the estimate 


u+ — Ur 


|o| 
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and thus 

In order to estimate this further, we need to determine the scaling of v. Using the 
estimate for the second derivative in ()in.5p with the fact that ReU has a maximum 
at Umax and no zero on the interval (^,Umax), we conclude that 

ReU(Umax)> 

Combining this inequality with the first inequality in ()11.27p . we find that the first 
term in (110.121) can be bounded in terms of the second term. More precisely, we obtain 
the inequality 

~ min(ei,e?) |U|2 ReA = C? jUp . 

Using this inequality in ()11.30p gives 

u (u+ — Ur)"^ ^ Gf . 

This concludes the proof on the interval (u, ^). 

It remains the consider the region (^, Umax)- We make use of the concavity of Re U, 
(llO.Sp . Denoting the zero of ReU by ui, we obtain 

ReU(Umax) ~ I Alp (ttmax “ Ulf . (11.31) 


On the other hand, for the potential to have a zero near the spectral parameter 
must scale like A ^ |Alp. Therefore, the second term in ()10.12l) can be estimated in 
terms of the first, so that 

i/3^min(e2,ei) |DpReU(Umax) = e2|UpReU(Umax). (11.32) 


Since < ReU(umax)/2 (cf. (|lU.12p and (I10.7p i. the estimate (|11.3ip can be extended 
to 


Re U( Umax) ~ |A1| (Umax ^) ^ ^ * 

Moreover, 

Re V'(u) ~ |Up (umax — u) for all u € [u^, u+] . 
We now combine the estimates (|11.32p . (|11.33p and (|11.34p to obtain 

V 


li-j- — U'p 


ReU' 


z^(u+ 


_ 1/3 _ e? |Up ReU(Umax) 

Re^ U' Re^ U' 

_ Gj |Up (Umax - ^ 2 

|A1P (Umax - It)^ 


(11.33) 

(11.34) 


This concludes the proof. 


□ 


We now estimate the Riccati solution in the Airy region. 

Proposition 11.11. Assume that in in the Airy region [ur-,u+], one of the following 
two conditions hold: 

(a) The potential has a positive imaginary part, ImU|[^J^_^J_^] > 0. 

(b) The imaginary part of the potential is small in the sense that 

I ImU] < |Up“^ for a suitable constant (5 > 0 . 
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Moreover, assume that the Riecati solution begins in the upper half plane, lmy{ur) > 0. 
Then there is a constant C2 (depending on Qi) such that for large |0|, the solution on 
the interval [ur,U-^-] can be estimated in terms of y{ur) by 


\y{u)\ < 62 \y{ur)\ 

lmy{ur) 

lmy{u) > --- 

*-2 


\4>{Ur)\ 

62 


< \(l){u)\ < 62 \(p{Ur)\ ■ 


(11.35) 

(11.36) 

(11.37) 


Proof. As in the proof of Proposition 111.91 we use the T-method choosing a and j3 as 
in (lll.24p . Then V and U are again estimated by (lll.25p . 

In case (a), we apply [121 Theorem 3.3] for <7 = 0. The error terms Pli,..., i?4 are 
estimated as follows, 


E 2 — P/4 — 0 

|^i|<^|ReP-ReP| + ^< 

/v V V 


l^sl 


^ |ImR| ^ |11| 


J ^ \/P ('U-l- — Ur) + 


Re V' + -Li (tt^ — Ur) 


/V 

= y/u (tt+ — Ur) H ^ Re R(u+) — Re R(ur)^ + ^ {u+ — Ur) 

<61 + 2 + 61 


As a consequence, we can apply Lemma lll.lOl to obtain 

1 /■“+ 

V 

1 
V 

N _ 

V 

Finally, by combining (|10.12p with (IIP.ID . (|10.2I) and (|10.7I) . we conclude that 1/ > |n|. 
This concludes the proof in case (a). 

In the remaining case (b), we choose g as 

g = . 

Then the error term PI4 is estimated by 

/*Pimax 


\EA\=g [ 

J u 


- I ImRl < ^ |^|i_, (n. - n.) CfHl ^ ^ ^ 1^,1-. 


VW\ 


a/CT 


This concludes the proof. 


□ 


In the next lemma we compare the imaginary part of the potential on the two Airy 
regions [u(i,u^] and 

Lemma 11.12. In the Airy case, one of the following three statements holds: 

(i) and > 0 
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Proof. We first consider the case that the Airy regions are near the poles in the sense 
that ri+,(vr — u^) < Then the factor sin^ u in (j2.8p is bounded by 

implying that Im(y + A) is bounded by Therefore, depending on the value 

of ImA, we are in one of the above cases (i)-(iii). 

It remains to consider the case that the Airy regions are away from the poles. Then 
the factors l/(sin^n) in (|2.8p is bounded by As a consequence, the imaginary 

part of V can be related to its real part by 

2 Tm Q 

ImV = Rel/ + (const) + 0(-\/|f^) . 

The function ReR has a zero in each of the intervals and It 

follows that at these zeros, the imaginary part of the potential has the form 


lniV{vL/R) = (const) + O(vl^) • 


It remains to show that on each of the intervals and the total 

variation of the function ImR is < By symmetry, it suffices to consider the 

interval [u^,u^]. For ease in notation, we omit then index L. We then obtain 


ImR < (u^ - u^) sup | ImR' 
[upu^\ 




|0| ™ |0| ^ y^ 

^ Cl 


This concludes the proof. 


□ 


11.6. Estimates on the Interval [u+,Uniax]- R remains to estimate the solution in 
the Airy case on the interval [u+, u^ax] (see (I10.15h and (llO.lOp ). 


Lemma 11.13. For any e > 0, hy increasing Ci one can arrange that 


Moreover, 


PU, 

J 11 + 


Wmax 2 

w-"' 


l^>WI 


> 




(ReR(M))4 


e “+ 


(11.38) 


Proof. In Proposition 111.11 it was shown that the WKB conditions (|ll.ip are satisfied 
on the interval [«+, Umax] ■ Thus the solution is well-approximated by the WKB solution 


<l> 


1 

(ReR)3 


Cl ^ + (72 e' 



(11.39) 


with error terms which are under control in view of the estimates in m- Note that 
one of the fundamental solutions in (jll.39p is exponentially increasing, whereas the 
other is exponentially decaying. 

Combining the estimate of Proposition 111.11 at u = Ur with the estimates (I11.35P 
and ()11.36l) on the interval [ur, u+] and taking into account the normalization (16.711 . one 
sees that the coefficient of the exponentially increasing fundamental solution in (I11.39P 
is bounded away from zero, and that \(j){u+)\ ~ |ReR(u+)|“4. This gives (lll.38p . 
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Next, we increase 65 to a new constant S 3 . Denoting the corresponding boundary of 
the Airy region by 12 +, we obtain 




> 


As a consequence, 

climax 2 

A+ W ^ 


(ReV) 

- 2 /“+Re^A^ 


1 r+ReVV f- ReVV 

-_ gJa+ 


: PUn 

J 14+ 




-2r+ReVvf^ -2r™Reyv\ ^ -2f“+ReVF 
= e ■'“+ 1 - e ■'“+ < e 


(11.40) 


The integral in the last exponent can be estimated from above by yT'(i2+ — u+). 
Applying Lemma 111.101 this term can be made arbitrarily large by increasing Si 
(and consequently S 5 ). As a consequence, the last exponent in (lll.40p can be made 
arbitrarily small. This gives the result. □ 


12. Integral Estimates of ImE 


In this section we shall derive the following estimates. 


Proposition 12.1. For all eigenvalues the following inequality holds, 

r" |ImE| ^ ^ 

Jue ^ ^ 


Proposition 12.2. For any d > 0, hy increasing Si one can arrange that that for all 
eigenvalues X in the WKB case and the parabolic cylinder case the following inequality 
holds: 



ImV 


< d. 


( 12 . 1 ) 


12.1. Elementary Estimates of the Potential. We begin with integral estimates 
of our potential. 


Lemma 12.3. The function ReE is monotone increasing on the interval [tto,Ur]- 
Moreover, 

nUr I I 

/ i~ • (^ 2 . 2 ) 

Juo \ReV \2 

Proof In order to prove the monotonicity of ReR, we first recall that in the proof 
of Lemma 110.11 we already showed that Re V is monotone increasing on the inter¬ 
val (u,Umax)- On the remaining interval [uo,u\, we need to consider the cases k = s 
and k ^ s separately. In the case k = s, the monotonicity of ReR is obvious 
from (lll.lip . In the case k ^ s, we see from (111.181) that ReR is convex. Com¬ 
bining this with the fact that uq is chosen as a minimum of ReR (see (111.171) 1. we 
conclude again that ReR is monotone increasing on the interval [uo,u]. 

For the integral estimate (112.21) . we first consider the interval {u£,U 2 ) with 

U2 = min (^Ur, 


(12.3) 
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Then the desired estimate is obtained by using that the integration range scales 
like . On the remaining interval (u 2 ,Ur) we consider the regions (u 2 , f) and (f,Ur} 
separately. In the first region, we approximate Re V by the quadratic potential 

I ReR| > a — c{u — ^ 2 )^ j (12.4) 


with parameters a, c > 0 to be specified below. Applying Lemma 19.11 it follows by 
explicit computation that 




' U2 


W < 

|R|^ 


J U2 


min(ur,f) 


(a - c{u- U 2 Y) 


u- U2 


■\/a-c{u- U 2 )" 


< 


U-U2 


miriK.f) « vClPl 


min('Ur,f) 


where in the last step we used that we are in the WKB region (see ()in.l5p . (jl0.14|) 
and (llU. 12 p i. Since the quadratic polynomial a — c{u — ^ 2 )^ is positive, we know that 


U — U2 < 



implying that 


L 


min(ur,f) 


1 


A < 

|R|i 


1 


y^CiOC |fl| 

Using that a > Ci |n| and c ^ |np, we obtain the desired estimate. 

On the remaining interval we estimate ReU on the interval [uo,Ur] by the 

quadratic polynomial 


I Re U| > a + b {u^ — u) + c {ur — tt)^ , (12.5) 

where the values of the positive coefficients a, b and c will be estimated below. It 
follows by explicit computation that 



< 



(^a + b (ur 


u) + C {Ur 
u) + C {Ur 



2 

y/ab + 2ay/c 


It remains to analyze the coefficients a, b and c. In view of (I10.12p . (jl0.14p and (jl0.12p . 
we can choose a > C 5 |n|. Moreover, the expansions ()10.3p - (|10.5l) show that we can 
choose either 5~|Oporc~|np. This concludes the proof. □ 


Lemma 12.4. For any 5 > 0 there is a constant C = C(5) such that the following 
inequality holds, 



1 C 

~W\' 


( 12 . 6 ) 


Moreover, there is a constant C such that 


r < c 

Jn, M 


(12.7) 
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Proof. As in the previous lemma, we first consider the interval {ug, U 2 ) with U 2 accord¬ 
ing to (I12.3p . Then, according to (|10.7I) and (|10.15l) . (110.141) . (|10.12ll and (Ill.Sp . we 


know that 


and thus 

n 1 ^ U 2 -U 0 

Juo ^ a/Ci 


The desired estimate is obtained by using that the integration range scales like . 

It remains to consider the interval {u 2 ,Ur). We consider the regions (u 2 ,|) and 
{^,Ur) separately. In the first region, the first and second derivatives of ReF are 
non-negative. We again estimate ReR by the quadratic polynomial (112. 4p . where the 
values of the positive coefficients a and c will be estimated below. It follows by explicit 
computation that 


rm 

J U2 


min(ixr ,J) 




< 


J U2 


min(nr,f) 


= —= arctan 
VC 


— c{u — U 2 Y 
^/c{u- U 2 ) 
■s/a - c{u- U 2 Y 


min(wr,5^) 


Since the arctan is bounded, it suffices to note that c ~ to obtain the desired 
1/1 n| behavior. 

On the remaining interval {^,Ur), we again estimate ReR by the quadratic poly¬ 
nomial (jl2.5D . where the values of the positive coefficients a, b and c will be estimated 
below. An explicit computation yields 


rur 2 P 

J U2 -x/lRl J U 2 


1 

7s 


U 2 y/a + b{Ur — u) + c{Ur - u/ 

= log ^6 -|- 2c {Ur — u) -|- 2^/c ^/a + b {Ur — u) + C {Ur — 


U2 


( 12 . 8 ) 


It remains to analyze the argument of the logarithm. From the explicit form of the 
potential (12.8h we know that we may choose 6, c ^ |np. If a ^ |np is large, the 
logarithm in (|12.8I) is given approximately by log(2-y/c ^/a). As a consequence, the 
difference of the logarithms at the upper and lower boundary points is uniformly 
bounded. Using that c |np, we obtain the desired estimate. 

It remains to consider the case that a < Then the arguments of the logarithm 
scale like |np, both at the upper and lower boundary point. As a consequence, the 
logarithm is again uniformly bounded, giving the desired estimate. This concludes the 
proof of the inequality (|12.6I) . 

In order to prove (112.71) , we first note that the estimate (112.61) fails to hold in general 
if (5 = 0. The problem is that in this case, the parameter b in (|12.8I) could be small, 
leading to a factor log |n|. This gives (112.7p . □ 


Applying Lemma [9T] to (112.71) . we obtain the estimate 
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Unfortunately, the factor log |n| is not good enough for our purposes. The next lemma 
shows that, if ImU vanishes at the right boundary point, then we get an estimate 
without such a logarithmic factor. 


Lemma 12.5. Assume that u G {u£,Ur) is a point where lm.V{u) = 0. Then 

r \imv\ 

Ju, ^ ^ 

Proof. Combining Lemma l9.II with Lemma 112.41 it remains to consider the case Ur > 
^ — (5. Moreover, it remains to estimate the integral over the interval We again 

estimate Re V by the quadratic polynomial (|12.5p with positive coefficients a, b and c. 
Moreover, we estimate ImU by 


We thus obtain 
ImU| 

' U2 

< 


r \ imv\ ^ r 

Ju, ^ Ju. 


I ImU(u)| < |n| {u-u) . 
|fl| {u — u) 


Q. r- 


^/a + b {Ur — u) + C (Ur — u)^ 
2c {Ur — u) +b n 


< 


pu, 

J U2 


|n| (Ur — u) 


^ya + b {Ur — u) + C {Ur — u)^ 

U2 


-—---^-= = —\JaPb{ur -u) + c(ur - u)^ 

y a + 0 (Ur — u) -h C (Ur — U)^ C 


The result follows because c |np and the argument of the square root is bounded 

by<|^^p. □ 


12.2. Estimates on the Interval (0,ui). 

Lemma 12.6. For any 6 > 0, we can arrange by increasing C 5 that 

rU£ 

Jo 

uniformly in and X. 


Proof. We begin with the case k = s. Using the asymptotics as in the proof of 
Lemma 111. 71 we find that on the interval (O,^^), 


\4>l{u)\ < 2|c| log(v^u) 

l&(«)l < T (aictan {loglv/fiu)) + |) < 


uniformly in /r. As a consequence, 

\4p\‘^ uniformly in /x . 

Applying Lemma l9.11 we obtain 






M M ej 

Re A C 5 


This gives the result. 
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In the case k ^ s, we work similarly with the asymptotics (111.221) . 

\4’\ ~ 


f 


1 


L „.2L 




/ |Imy||<?ig|2<|0||H^nf+2 < 

JO 


mo. St g 2 L +2 ifTo^ 


ReA| 


This concludes the proof. 


□ 


12.3. WKB Representation of We again let (f)^ be the solutions (|2.17p with (j)L 
defined by ()10.17l) and (|10.16p ). In this section, we shall approximate (pY, on the 
interval {ui,Ur) by a suitable WKB wave function. Our starting point is the WKB 
approximation of 4 >l 

^TTTTy-T^ r’^ . /T7 

(12.9) 


, ^ , CWKB r ^/v 

(pL ~ 0WKB = —e-'“o 


In order to comply with the initial conditions (16.7p . we must choose 




"WKB 


Imy^K(Mo) 

In the next proposition, we compute what this approximation means for the osculating 
circles as introduced in Section [6l 

Proposition 12.7. On the interval [ui,Ur], the radius and center of the oseulating 
cirele are given by 


P ~ C('wo) + 


1 


2 c 2 

^OwKB 


and 


„ 1 -2 r ReVV 

Rk, - e ■'“0 
2 


( 12 . 10 ) 


where the error can be made arbitrarily small by increasing Ci. 
Proof. Using the WKB approximation (|12.9p . we obtain 

U' 


y = 


Au- 


4U 


and thus 


({u)-C{uo)=[ \/U(r) e 


1 


2 c^ 


f 

J un 


d -2 r w , 1 

-C JUQ — - 


dr 


2 c 2 

^'-WKB 


-2 vW ^ 
e ■'“0 — 1 


( 12 . 11 ) 


■-WKB JUQ 

It is remarkable that the integral can be carried out explicitly, giving a simple expres¬ 
sion for C{u). □ 

We next compute (fY- Using (|12.9I) and (112.lip , we obtain 

A) = (t>L{u) Cl{u) = ^ + (C(A - C(«o))^ 


^WKB UT 

</V 


Jo ^ 


2 c 2 

^^WKB 


-2 r W . 

g JuQ _ I 
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Hence 




2c 


WKB 


a r w 1 - r w 

_ pJUQ * _ _ p '^'^0 


( 12 . 12 ) 


where a is the constant 


CX — 1 “h 2 c 


'WKB 


fUQ 

Jo 


Lemma 12.8. By choosing 65 sufficiently large, we can make the expression 

\a\ — 1 

arbitrarily small. 

Proof. Choosing 65 sufficiently large, we can arrange that the potential at uq is ap¬ 
proximately real and negative (cf. (|11.12l) . (111.141) and ()11.19p . (111.1811 ). Hence 


^WKB ~ * ) 


implying that 


a 


1 + 


r^o i 

^ 7 o ^ ^ 2 iC(uo)| , 


(12.13) 


with an arbitrarily small error. Next, using the initial conditions (|6.7I) in (16.101) . we 
obtain 

% 1 
P{uo) = C{uo) - - and R{uo) = - . 


2 ^ 2 
Solving for ({uq) and substituting into (112.13p . one finds that |q!| 


2 |p(uo)|, and thus 


a| — 1 


< 2 ||p(no)| - i?(uo)| . 


Since C(0) = 0, we know that |p(0)| = i?(0). Hence 

ruQ 

\\p{uo)\ - R{uo)\ < (Ip'I- h |i?'|) . 

Jo 

In view of (I6.10p and (|6.12l) . we know that |p'| = Hence our task is to show that 
the total variation of R on the interval (0, uq) is arbitrarily small. Since R{uq) = it 
suffices to show that the total variation of logi? is small. Thus, according to (j6.12l) . it 
remains to estimate the integral 


ruo |^/| ^ 

Jo ^ Jo 


Imy 


In view of the estimates in Sections 111.21 and 111.31 we know that 

Im y > Re . 

Hence 


/ 


|ImH| < < \n\ 


Im y Re /X Re /r 
By choosing C 5 sufficiently large, we can make this expression arbitrarily small. 


□ 


Using this lemma in formula (jl2.12p . we obtain the estimate 




< 




cosh 


2 f ReVV) . 

JUQ J 


(12.14) 
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12.4. Integral Estimates of WKB Solutions. 

Lemma 12.9. Assume that u E {u£,Ur) is a point where lm.V{u) = 0. Then 



< 1 . 


Proof. We can work with the WKB approximation (112.141) . According to Lemma [12.5l 
we know that 



< 1 


giving uniform control of the absolute value of the hyperbolic cosine in (112.1411 . As a 
consequence, 


ru 

J U£ 


D|2 




pu 

Jui 


\lmV\ 


<1 


where in the last step we again applied Lemma 112.51 This concludes the proof. □ 


Next, we take the absolute square of the WKB approximation (112.12p . 

1 5 Rp^/V 1 1 


\4>Liu)[ 


\a\ 


4 I CwKB I ^ ■\/\V\ 

1 1 

+ 


2 r ReVV , 
e ■'“0 -|- 


4|cwKBp \/\V\ 


-2 r ReVV 
e '^'^0 


2 I CwKB I ^ \/\V\ 


T-, / 2i r“ Im vT 

Re a e ■'“o 


(12.15) 

(12.16) 


The integrand of the last term is oscillatory. As a consequence, the resulting integral 
is small, as quantified in the next lemma. 


Lemma 12.10. For any 5 > 0, by increasing Ci one can arrange that that for all 
eigenvalues A the following inequality holds: 


r- |ImK| 


2i ff Im W 

g ■'^0 


< <5. 


Proof. Since on the interval (0, ^), the function Im V changes signs only once, it suffices 
to show that for any ni,tt2 E [ui,Ur], 


rU2 

J U-[ 


Im V 2i Im y/V 
— :- e -^^0 


<5. 


Integrating by parts. 


I- 

J U ^ 


ImK 2 i r ImW 

— - p *^1^0 — 

vWJ 


J U ^ 


ImK 


d 2i P Im pV 


Y^|l/| 2i Imy/V du 


ImV 


2i f“ Im vT 

JUQ 


y/\V\ 2i Im\/V 


U2 

p d / 

Ui 

Jui du \ 


ImV 


du y y^|K| 2i \my/V^ 


2i r Im W 
e ■'“0 
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we obtain the estimate 


rH 
J U\ 


“2 IvnV 2 i/“ ImW 

- g •'“0 


< 


< 

rs_/ 


Imi/| 

2 1^1 

ImFl 

|ReF| 


+ 


Ul 


+ 


ImRl 

2 1^1 

ImRl 

|ReF| 


+ 


U2 


rU2 

J Ul 


du 


ImR 


2y/\V\lmVV 


+ 


r“ 2 |iniR'| r“2|ImR'|2 \ReV'\\lmV\ 

4 WW^Jui 


ReRp 


All the terms except for the last summand can immediately be estimated in the desired 
way using the explicit form of our potential. For the last term we use the monotonicity 
of Re R (see Lemma 112.31) to obtain 


J U^ 


ReR'l |ImR| 


| 0 | / 

J Ul 


I Re R'l 


|ReR|2 


,,, |ReR|2 
This gives the result. 

Keeping track of the constants, we now write (JP as 


< |n| sup |R| ^ 

[ui,U2\ 


(tP{u) = 


cl4>1 


if u < ^ 


CR (j)^ if u> 


□ 


(12.17) 


2 ’ 


where cr and cr are non-zero complex numbers (and and are again the solu¬ 
tions (I2.17P with (pL and (pR defined by (I10.17h and (llO.lbp ). 


Lemma 12.11. If\cL/cR\ < 1, then 

r \luiV\\cP^\^ <\cr\\ (12.18) 

■^2 

Proof. We again denote the zeros of ImR hy ur and ur. Lemma [l2.6l and Lemma [l2.9l 
imply that 

Moreover, using the representation (112.171) in (19.21) . we obtain 

ruR 

I \lmV\\cP^\^<\cR\\ 

2 

This gives the result. □ 


Lemma 12.12. If\cR/cR\ < 1, then 

r-' |ImR| ^ ^ 

Juo ^ 

Proof. The strategy is to combine Lemma 112.111 with the fact that the potential is 
approximately symmetric with respect to reflections at u = ^. In order to make this 
approximate symmetry precise, we consider the homotopy 

Rt-(u) := r R(u)-|-(1 — r) R( 7 r — u) for r G [0,1] . (12.19) 
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Then the mean value theorem implies that 

I Rel/i(tt) — Re V 2 (u)| < sup -^ReVr{u) , 

re [0,1] dT 

and similarly for the imaginary part. Using that the function sin^ u in (|2.8I) is reflection 
symmetric, one finds that 

|ReUi(u)-ReU 2 (M)| < and | ImUi(u) - ImU 2 (u)| < 1 . (12.20) 

This implies that the WKB approximation holds on the “reflected WKB-region” [vr — 
Ur,7r — uo]- Using the WKB approximation (112.15p and (112.lOp in (112.18p . the oscil¬ 
latory contribution (112.lOp was estimated in Lemma 112.101 Noting that one of the 
factors exp(±2 Re\/U) in (|12.15l] is greater than one, Lemma 112.111 implies that 


rn—UQ 
J IT —Ur 


|ImU| 


< 1 


Again applying the reflection argument and the mean value theorem, it remains to 
show that 

d /Iml4^ ^ ^ 


d nmVr\ 
Juo dr V VVr ) 


Again using the explicit form of the potential 


d /ImK\ 


IrndrVr 

1 

{drVr) ImVr 

dr V yAf ) 


y/Vr 


3 

u/ 


< 


1 ^ | 0|2 


VW\ |u|t ■ 


Integrating this inequality from uq to Ur, we can apply Lemma 112.31 to obtain the 
result. □ 


Proof of Provosition \12. 1[ It suffices to consider the case \cl/cr\ < 1 and to show that 



Then case {cl/crI > 1 can be treated similarly by exchanging the left and right 
subintervals with the reflection u tt — u. 

On the interval [u^,uf’], we argue as in the proof of Lemma ll2.11l to obtain 



On the interval [uj^,u^] we consider different subintervals: The region from ui to uq is 
estimated in Lemma 112.41 The region from uq to Ur, on the other hand, is estimated 
in Lemma ll2.12l This concludes the proof. □ 
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12.5. Estimates in the WKB and Parabolic Cylinder Cases. We now give the 
proof of Proposition 112.21 Our strategy is to refine the method of Section [9] and to 
combine it with the “reflection argument” which was already used in the proof of 
Lemma 112.121 Another ingredient is Proposition 112.11 (whose proof was completed in 
the previous section). 

We will apply Proposition 112.11 in the following way. In the WKB region, we know 
from (jlO.ip . (|10.2I) and (|10.14p that 


and thus 


ImF|<|L!| 
|ImK| 


<1 


and ReK<-ei|ll| 


in the WKB region . 


I Re K| ~ Cl 

As a consequence, we may expand the square root of the potential as 

ImV 


W = VReV + ilmV = i^\ReV\ -ilmV = i^/\ReV\ + 


vr^\ 


+ 


showing that 


In particular. 


Re^ = 


ImV 


Re^l < 


|lmK| 


( 12 . 21 ) 


’ 

and applying Proposition 112. we conclude that 


rUr 

/ |Re^|<l. 

Jui 


( 12 . 22 ) 


This shows that that the exponentials and hyperbolic cosine in the WKB approxima¬ 
tion (see (|12.9p . (|12.12h . (112.14p and ()12.15h . (I12.16h i are uniformly bounded. 

We again assume that A E C is an eigenvalue and the corresponding eigen¬ 
function. Moreover, assume that we are in the WKB case or the parabolic cylinder 
case (but not in the Airy case, which is excluded in Proposition 112.21) . Using the 
representation (112.17p in (19.2p . we obtain the identity 


|cl| 


/*7] 

"L 


ImU = 0 


ImU |</.°|2 + |cK|2 
JO 

Denoting the integrands by /, 

f{u) := 

we decompose / into its even and odd parts, 

/ = /+ + /- with f±{u) := i (^/(u) ± /(tt - n)) . 


(12.23) 


ImU \ 4>i\^ if u < ^ 
ImU if u > I , 


(12.24) 
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Then we can rewrite (|12.23l) as 


0 = \ cl \ 


= CL 


TT 

'' [\u 

Jo 

TT 

'' [\u 

Jo 


+ /-) + |c_r|^ / {f+ + /-) 


+ /-) + |c_r| 


(/+-/-) 


= ( IclP + |c_Rp 


/+ + ( |cl| - |ci?| 


/- 


Since the constants cl and cr are non-zero, we obtain the inequality 


u 


< 


1 

i: 


f- 


Combining this estimate with ()12.24l) . we obtain the inequality 


£ 




D|2 


TT 

[\u 

Jo 


+ /-) 


< 2 


£ 


f- 


(12.25) 


We now estimate the right side of this inequality obtain the following result. 


Lemma 12.13. For any 6 > 0, by increasing Ci one can arrange that that for all 
eigenvalues A in the WKB ease and the parabolic cylinder case the following inequality 
holds: 





< 5 . 


Proof. We introduce the “parity transformation” P which reflects at the point 

P : (0, tt) —>■ (0, tt) , Pu = TT — u . 


Then (112.25p can be written as 


Jo 


< 



Im(yoP)|())goP|") 


On the interval (0, this integral can be estimated by Lemma 112.61 where we 
choose 5 = (5/16. Moreover, in view of the estimates of Lemma 111.91 in the para¬ 
bolic cylinder case we may estimate the functions (fP and C, at the point Umax in terms 
of their values at Ur- This shows that the integral over [u^, §] can be made smaller 
than (5/8. We conclude that 


TT 

Jo 

We next specify the wave functions and (f)^. Using (jl2.15|] and (112.16|] together 
with Lemma 112.81 and Lemma 112.101 we know that in the integral on the right side 
of (112.26h . the factor |(/>°P may be replaced by the function 

—; - cosh \ 2 £ Re \/U ) , 

2v^ V Ao ) 

making an error which can be made arbitrarily small by increasing Ci (here we use 
Proposition 112. Il to conclude that a small pointwise error gives rise to a small error of 


< 


f VimU |(/.gp -lm(UoP) \ f^oP\ 
J UP ^ 


(5 

+ 4 


(12.26) 
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the integral). Moreover, using (jl2.2ip . we may replace KeVV by ImV/■\/\V\^ again 
making an arbitrarily small error. Therefore, setting 


P^{u) := 


2v1^ 


cosh 2 


r lmV\ 

Juo 7W\) ’ 


we can arrange that 


ru, 

J Ufi 


Im F I (/ij r — / Im F 


ru, 

J ue 


5 

< - . 
- 8 


(12.27) 


(12.28) 


Using the same argument on the interval [f,vr], we conclude that 


7 

i: 


< 


lmV\(j) 


D|2 

lI 


< 


L 


r (imV pl-lm{VoP) (pgoP)) 
Jui ^ ' 


+ 


Im(U- UoP) pg 

r\m{VoP) U-{p^^oP)) 

J Up ^ ^ 


+ 


(12.29) 


(12.30) 


We next estimate the integrals in (jl2.29p and (112.30p after each other. In order to 
estimate (112.29p . we first note that, from the explicit form of the potential (|2.8p . it is 
obvious that 

I Im (U - U o P) I < 1 . 

As a consequence, 

rUr PUr PUr 1 

/ \lm(V-VoP)\\p’l\< \pI\< 

Jug Jui Jui y \ v \ 

where in the last step we again used Proposition 112. Il to conclude that the hyperbolic 
cosine in (|12.27p is uniformly bounded. Using the estimate (112.71) in Lemma 112.41 we 
conclude that 

Im(U-UoP)||pg|<^°^l^l 


PU, 

Jui 


llll ’ 


which tends to zero for large |11| and can thus be made smaller than 5/A. 

In order to estimate (|12.30p , we again use the homotopy (112.191) . Setting 

uo(r) = tuq + (1 - t)uo 


Pt = 


1 


‘A^/Wr 


cosh 2 


PU 
J UQ 


ImVr 


(t) yWr 


we again use the mean value theorem to obtain 


\pl -{pr°P)\ = \pi - /5o I < sup 

Te[o,i] 


d , 


^ sup 


1 




\drVr\ |ImK 


|5rUo('r)| + f 

Juo{t) 


|(9t- Im U-1 ^ \dr Im U- drVr 


+ 


lK|i 


(where we again used (|12.22p to conclude that the hyperbolic cosine is uniformly 
bounded). Using that (see also (I12.20p i 

|5rReU| < |11| , |a^ImU|<l and |5^uo(r)| <-^ , 
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a straightforward computation shows that for any 5 > 0, we can arrange that 

5 


\Pl - (prO^)I < 


Hence 


r\m{VoP) te-(pSoP)) 

Jui ^ ^ 


< 6 


VW\' 

r- |Im(HoP)| 

Jue 


<<5, 


where in the last step we again applied Proposition 112.11 By choosing 5 sufficiently 
small, we can arrange that 


r\m{VoP) te-(pSoP)) 

Jui ^ ' 


6 

< - . 
- 4 


This concludes the proof. 


□ 


Proof of Proposition \12.^ The remaining task is to estimate the integral in (|12.1I) from 
above by the integral in the statement of Lemma 112.131 Applying (112.28p , we obtain 


ru, 

J ue 




D|2 


> 

rUr 

/ ImHpg 


J Ul 


-(5. 


Moreover, using (112.271) . we obtain 


rUr PUi 

/ lmVpl = 

J Ul U£ 

1 


ImV 


cosh I 2 


r ImP\ _ 1 r 

Juo ~^Ju 


d 


2v^ V" Ao VW\J 4 du 


= — sinh ( 2 
4 


f 


ImP\ 1 , , 

I-smh 2 


'no V\V\J 4 

In the last summand we can use the estimate 

\n\ 


^ / r ImH\ 

V Juo yW\J 

ImP\ 


ImP\ 

Juo VWj) ' 


rut 

J uo 


< 


~ yelM ~ v/ET 


to conclude that this summand can be made arbitrarily small. In the first summand, 
on the other hand, we apply the inequality | sinhx| > |x|. This gives the result. □ 


13. The A-Dependence of the Osculating Circles 

In view of the result of Proposition 18.11 it remains to consider the situation for 
large |n|. In this regime, Weyl’s asymptotics as worked out in Lemma 17.81 is of no use, 
because we have no control of how the error terms 0(n) and O(n^) depend on Q. In 
particular, we cannot expect that the gaps between the eigenvalues for a real potential 
are so large that the imaginary part of the potential can be treated as a slightly 
non-selfadjoint perturbation. As a consequence, we must analyze the spectrum of 
the Hamiltonian with the complex potential (j2.3l) in detail. As a technical tool, we 
will again work with the osculating circle estimates as developed in Section O As 
a refinement, we need to analyze in detail how the osculating circles depend on the 
spectral parameter A. In view of (12.8p and (j2.9p . we know that 

Ha := 5aH = -1. (13.1) 

Moreover, we choose uq in ()6.8p independent of A (where yo clearly depends on A). 
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In the following computations, we treat Re A and Im A as two independent real 
variables. Then for any function /(A) which is complex differentiable, we have 

r f = wv—r / = hm--- = ifx ■ 

a Re A a Im A h\o h 

In the next lemma we compute the A-derivatives of 4>. 

Lemma 13.1. Choosing the initial conditions (16.71) . 

d , ^ Irayxiuo) , 

log (p{u) = - -^^ + 


9ReA 

d 

5Im A 


2Imy(uo) 


/ 

J U, 


y\ 


Reyxiuo) . 


no 


L 


y\ ■ 


no 


Proof. Differentiating (16. 7p and (I2.15p with respect to A gives 


d 


dReX 

d 

5Im A 
d d 


logfiuo) = - 


1 Imyxjuo) 

2 Imy(uo) 


1 X 1 Reyxiuo) 

logcpiuo) = -- 


log fiu) = 


2 Imy(uo) 

yiu) = y\iu) 


dudRoX 9 Re A ‘ 

Integrating the last two differential equations from uq to u gives the result. 
We next compute the second mixed derivatives of K and p. 

Lemma 13.2. Choosing the initial conditions (EH), 
d d , ,, ImR Imy>, 


□ 


(9 Re A (9u 
d 


\ogK{u,X) = -- 
d 


dReXdu 
d d 


piu,X) = -i 


Im^ y 

ImV flmyx Imyxiuo) 


d Im A du 
d 


log K{u,X) = - 


(fPlrr? y \lTa.y 2Imy('Uo) 
1 ImR Reyx 


+ 


nu 

/ yx 

J no 


Imy 


Im^ y 


d 


d Im A du 


piu,X) = - 


2 (fp Im^ y 


— i 


ImV (Royx ReyA(ito) 


(fPlviT y \lTiiy 2Imy(uo) 


nu 

+ i yx 

J no 


(13.2) 

(13.3) 

(13.4) 

(13.5) 


Proof. Comparing ()6.13l) and (|6.14l) . one sees that 

— logR:(u, A) = -- 

ou Ray 


Differentiating with respect to A gives 

d d 


dReXdu 
d d 
d Im A du 


logK[u,X) = 
logKiu,X) = 


d Im V 
9 Re A Imy 
d Im R 
film A Imy ’ 
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and using (113. ip gives (113.21) and (113.411 . In order to derive (|13.3p and (|13.5I) . we first 
apply (16.101) to obtain 


K 


d d d 

p{u,X) = K 


d'R.eXdu 


K- 


d 


= le 
d 


-2i'd 


ImV 

Imy 


i ImV 


(9Re A lur y 

^9 „ImyA 

-2 „„ , log(/> - 2 -- 

imy 


d Im A du 


p{u,X) = K 


9 Re A 
d f i IraV 


= ze 




1 


Imy 


(9 Im A Ivc? y 

-l-2Iml/ ^ 


9Im A 


log cj) — 2 Im V 


Re 1/A 
Imy 


Applying Lemma 1 13.1 1 gives 
d 


K 


K 


9 Re A c/m 
d d 




ImyA(Mo) r 
2 Imy (mo) Juo^^ 


p{u,X) = - 


Imy \ Imy 

^ ^KeyA ReyA(Mo) 


— 2ie 


d Im Xdu^'' ^ ' Im y 
Using ()6.10l) and (|6.5I) gives the result. 


Imy \ Imy 


2Imy(Mo) luo^^ 


□ 


13.1. General Estimates in the WKB Region. In the previous section, we derived 
general formulas for the A-derivatives of the center and radius of the osculating circles. 
We want to use these formulas in order to control the behavior of the osculating circles 
in the WKB region (m^, m^). We now work out the corresponding estimates in general, 
stating the assumptions needed for the estimates to work (see eqs (113.6h . (113.7h . (113.8h 
and (jl3.9p below). We also derive explicit formulas for the error terms (see eqs (I13.1ip ~ 
()13.13p below). The remaining task will be to justify the above assumptions and to 
control the error terms. This will be done subsequently in the different cases in the 
following Sections 113.21113.31 

Consider an interval I = [ ua , Ub ] C (0,^]. We assume that on I the following 
inequalities hold for suitable constants C, c > 1, 



ImU| 

Imy 


< c 


< Imy < C a/I ReU| . 


(13.6) 

(13.7) 


Lemma 13.3. Assume that uq & I and that (113.6p and (113.71) hold. Then the func¬ 
tion K (m) (cf. ()6.10p ) is bounded on the interval I by 

- < K{u) < c , 
c 

where c = e^. 


Proof. Comparing (16.101) with (16.7p . we know that K{uo) 
and ()6.1ip we obtain 

d , K' ImV 

^ log - X “ • 

Integrating and using (|13.6p gives the result. 


2. Moreover, from ()6.10l) 


□ 
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We now estimate the terms appearing in Lemma 113.21 We again consider an inter¬ 
val I = [ua,Ub]. We assume that on I the inequalities (113.61) and (|13.7I) hold and that 
the following inequalities hold for suitable constants C,c> 1, 



|Rey| ^ c 

V\^iv\ - V|ReR(u)| 

ReR(u)| < inf |ReR| . 

[UaM 


(13.8) 

(13.9) 


Lemma 13.4. Under the assumptions (|13.6P " (jl3.9p . there is a constant ci = ci(c, C) 
such that on I, 

More precisely, 


y\{u) - 


2Imy(u) 


< Cl 


+ I^Al/ol 


+ Cl \/|ReR(u)| f 

Ua I 

Proof. We integrate by parts to obtain 


/|ImR| 
V\ \ Imy 


+ I Rey| 


(13.10) 



Using this relation in (17.131) gives 


y\{u) 


4)^{ua) f^ 1 \ 

(tP{u) V 2i lm.y{ua)) 

1 2i r ImU 

_/ —— Re 1 /_ 

2iImy(M) Ju^lmy \ 4Imy 


We now estimate the resulting terms: 




dxyo 


(jP^jUg) 1 
(/)2(u) 2i Im y{ua) 


1 

2y{u) 


. Imy(u) 1 /^ 2 - 
< c T-^—r I^Ayol <C c 


< c 


lmy{ua) 
lmy{u) 

Im^yiua) ~ |ReVb| 


|Qa?/o| 

vWWl 


a/| ReU(u)| < C^c \dxyo\ 


< 


ch 


A/|ReU(u)| < 


C^i 


A/|ReU(n)| 


/ C 1 
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Moreover, 


2 

r / 

(jp{u) 

Jua Imy \ 


Rey - 


< 2C^c^ 


\/|ReR(n)| [ 

J U 


ImR 
4Imy 
1 


< 2c^ Im y 


r 1 

Jun. 


Rey - 


ImV 


4Imy 


|ReR| 


Rey - 
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4Imy 


^ 5 ^^3.2 A/|ReR(n)| ^ 5 cC^c^ 

- 2 ^ ^ inf/ I Re R| “2 ^|ReR(M)| 

Combining all the terms gives the results. 


□ 


Proposition 13.5. Under the assumptions (jl3.6p ~ (ll3.9h . the following inequalities 
hold, 


d 


5Im A 


r 

logK(u.\)\l + 

J u 


^6 


d 


dReX 


log K{u,X) 


Imy 

lib 


dp{u,X) 

Ub 

L 

dp{u,X) 

Ub 

dReX 

'^a 


9Im A 

R-a 


where the error term £ is given by 
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(13.13) 


with a constant C 2 = C 2 (c, C). 


Proof. We integrate the formulas of Lemma 113.21 from Ua to Ub- Possibly after inte¬ 
grating by parts, we apply ()13.7p - (|13.9l) and Lemma 0.3.41 More precisely, we estimate 
the individual terms as follows: 
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Combining all the terms gives the result. 


ReR| |ReR| 

+ I^Ayol 


□ 


13.2. Estimates on the Interval [uo,Ur]. We shall now apply the estimates of Sec¬ 
tion [T3R] on the interval [uQ,Ur\ (see (|10.14l) . ( 111 . 121 ) and (|11.17l) i. Our task is to 
show that the inequalities (I13.6p . (113. 7p . (113.81) and (|13.9I) hold. Moreover, we need to 
estimate the error terms £ in (I13.1ip - (I13.13I) . We begin by collecting a few properties 
of our potential. 

Lemma 13.6. At the points uq and u (see (111.121) . (|11.19l) and (IIP.lip ). for suffi¬ 
ciently large 65 the potential satisfies the inequalities 

|ReR(uo)-ReA| < |11| (13.14) 

|R'(no)|<|0|i, |R"(uo)| < (13.15) 

ReR(no) < ^ ReR(u) < 0 . (13.16) 

8 

Proof. The inequalities ()13.14p and (113.ISp follow immediately from ()11.12p . (Ill.lip 
and (I11.19p . pil.lSp . Combining (113.14p with ()10.2p . we find that for sufficiently 
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large 65 , 


Rey(Mo) < Re A. 


Moreover, similar as in the proof of Lemma llO.il again for large 65 one obtains 


16 


ReR(li) > tt; ^ - ReA = - 


16 Re A 4c|0|^ 


17 


13 Re A 
17 


17 4 

^#>-hReA. 

Re A 17 


Re A 


— Re A 


Combining these inequalities gives (jl3.16h . 


□ 


Lemma 13.7. Assume that the condition (jl3.6p holds. Then the conditions (I13.7p . 
(I13.8|] and (jlO.Op are satisfied on the interval [uQ,Ur\. Moreover, choosing ci and C 2 
sufficiently large, we can arrange that the error terms £ (113.111) - (113.131) are bounded 
by 

£ < — / , . (13.17) 

20 Ju, vIR^ 


Proof. The inequality (113.91) follows immediately from Lemma 112.31 Since the WKB 
conditions are satisfied by Proposition lll.il we know that 

(13.18) 

The results of the analysis in mm gives us rigorous bounds for the error of this 
approximation. This implies (jl3.7l) for sufficiently large C. 

In preparation for proving the other inequalities, we need a few estimates for the 
potential. First, from the explicit form of the potential (12. 8 p . it is obvious that 

V'l IT/'"I < lr>|2 


W'UV" 


and |ImP'|<|0|. (13.19) 

Estimating the potential without derivatives is a bit more subtle because the constant A 
comes into play. We first note that by construction of Ur (see ()10.14p . (I10.12p . ()10.7I) 
and () 10 . 2 p L we know that 

ReR<-ei|f2|. (13.20) 

Moreover, using (|10.ip in ()9.3p . the imaginary part of the potential is uniformly 
bounded, 

|ImP|<|L!|. (13.21) 

Therefore, by increasing Ci we can arrange that the real part of the potential dominates 
its imaginary part in the sense that 

|ReR| 


ImE| 


> Cl . 

rs_/ ^ J- 


(13.22) 


We now come to the proof of (|13.8p . We first need to estimate the real part of y. 
Using (|13.22l) . we may express (|13.18l) in terms of the real and imaginary parts of V 
to obtain 

;-^- V j - Im V V 
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(where the errors are again under control in view of (113.22p and the fact that the WKB 
conditions are satisfied). The last estimate implies that 


r iReyl 

.0 


r\imv\ r \v'\ 

Ao |ReR| Juo\V\y 


(13.24) 


In the first integral on the right, we apply (|13.9I) to obtain 


r |ImI/| ^ C r \lmV\ ^ 1 

Juo iReRl - y^|Rey(?r)| Ao ~ v'|ReR(n)| ’ 

where in the last step we again used the WKB approximation (|13.18l) together with 
the inequality (113.61) . The second integral in (113.24p can be estimated by 

r \v'\ ^ r |K'| ^ r ReW r iimWi 

Juo |K|i ^ Juo |ReK|i ~ Juo |ReK|i Juo |ReK|i 

r/ 1 y \n\ ^ 1 | 0 | 

~ Juo ^ |ReK(u)|i ~ v'|ReK(u)| ^ \n\ y|ReK(u)( ’ 

where in the last line we also used again the monotonicity statement of Lemma 112.31 
together with (|13.20p . Combining the obtained inequalities gives pi3.8p . 

In order to prove (|13.17p . we first estimate the term dxyo. Differentiating (I11.13P 
and using the monotonicity of Re V as well as the results of Lemma 113.61 we obtain 


I^Ayol < 


1 

V\^eV{u)\ ' 


For this reason, we may disregard the term dxUo in ()13.11l) and (|13.13p . In order to 
estimate the integral (113.lip , we simply bound the first bracket in the integrand by 


Iml/| 1 

ReK| ~ ^ ■ 


(13.25) 


The first summand in ()13.12l) can be the estimated with the help the monotonicity 
of Re V by 


1^9 8 9 9 

|ReK(u^)| - |ReK(u,.)| “ |ReK(u)| “ |ReK(u^)| “ |ReK(uo)| 

rur I y r- ReV ^ r- 10 \V'\ 

~^Juo “Vno R^-Ao 

where in the last step we used again that the real part of V dominates the imaginary 

3 

part and that ReK is monotone. Now the factor |K^|/|K |2 is small in view of the 
WKB property established in ProDosition lll.il 

To estimate the integrand in (jl3.12l) . we fist use (113.231) . 

|Rey| ^ |ImK| |K'| ^ 2 r|ImK| |K'| ' 

|ReK| ~ |ReK|l ^ |ReK||K| ^ ^ ' 

Now the square bracket can be made arbitrarily small by applying again (113.25p and 
Proposition 111.11 
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It remains to estimate the nested integral (I13.13p . We first exchange the orders of 
integration, 

Im^V I Im I/'l + I Re y| I Im F| 


rUr 

/ du 
J UQ 
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dr 
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I Im V'\ + I Re yl I Im V\ 
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du . (13.26) 


ReR|5 

Now we estimate the inner integral term by term. The first term can be estimated 
with the help of (113.211) and Lemma 112.31 

1^2 y 
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Jt |ReR|i ^ Jr 
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ReR|2 ~ Cl 


Next, using (|13.19l) and (|13.20p . we find that 
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'ImR'l ^ 1 , , ^ 1 

< — [Ur- t)< 


ReR| ~ Cl 


Cl 


In order to estimate the remaining term involving Re y in (113.260 , we first apply (113.231) . 
Then the first summand on the right of (113.23p gives rise to a contribution which is 
precisely of the form of the first summand in the curly brackets in (113.131) . Hence it 
remains to consider the term 


i: 


\lmV\ \V'\ 


ReR| |R1 
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ImR| 


ImR'l +ReR' 


where in the last step we applied Lemma 112.31 The first summand on the right can 
be estimated with the help of (113.2ip and (|13.19l) by 
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where in the last step we again applied Lemma [12.3l The remaining second summand 
on the right is estimated as follows, 


r 






ReR' 
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|L!| 


<1 


|R|2 (ReR)2 - |ReR(nr.)| ~ Ci 

where we again applied Lemma 112.31 This concludes the proof. 


□ 


13.3. Estimates on the Interval [u£,uo]. In preparation, we note that, using (110.11) 
in (|9.3I) . we again obtain the following uniform bound for ImR, 

|ImR|<|L!|. (13.27) 

We treat the cases k = s and A: / s separately. 

Lemma 13.8. (Case k = s) Assume that k = s. On the interval \u£,uq], the con¬ 
ditions (|13.6p - (|13.9p are satisfied. Moreover, by choosing Ci sufficiently large, we can 
arrange that the error terms £ (|13.1ip - (|13.13l) are bounded by 
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Proof. Near the pole, the potential has the following asymptotic expansion (cf. (I11.12P 
and ()11.14l) l. 


Re 1/ < —- |n| — Re/i + 0(1) < — — Re^u < — Re/x < Re A (13.28) 

V'{u) = ^ + + 0(|fl|M) + 0(|fl|\^) . (13.29) 

Since (|13.29p has a positive real part, we know that ReR is monotone increasing. This 
implies (|13.9I) if we also keep in mind that the imaginary part of V is dominated by 
the real part in view of (|13.28p and (113.191) . Moreover, using the results of the analysis 
in mm, we know that (113.18p holds with rigorous error bounds. This implies (|13.7I) 
for sufficiently large C. The inequality (|13.6p follows from the estimate 
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where in the last step we used (jl0.2p and (jll.l2h . In order to prove (113.811 . we again 
apply (113.241) and estimate the two resulting integrals by 
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In order to estimate £, we proceed exactly as in the proof of Lemma [13.7l using the 
following estimates: 
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This concludes the proof. 
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□ 


In the case k s, the function Re R is monotone decreasing on the interval [u£, uq]. 
Therefore, when applying the estimates in Section 113.II we need to proceed backwards 
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in u, starting from uq. Therefore, the conditions (I13.8P and (113.91) need to be replaced 
by 


iReyl ^ c 

Ju '7W\~7WW\ 

ReF(u)| < inf |ReR| . 
[w,«o] 


(13.31) 

(13.32) 


Lemma 13.9. (Case k ^ s) Assume that k ^ s. Then on the interval [uh^uq], the 
conditions (113.6p . (113.7p and (|13.31l) . (|13.32l) are satisfied. Moreover, by choosing Ci 
sufficiently large, we can arrange that the error terms £ (113.Ill) " (|13.131) are bounded 
by 
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1 r° 1 

“ ^ Ju, ^ ■ 


Proof. Near the pole, the potential has the following asymptotic expansion (cf. (111.201) 
and (jll.lSp h 


ReR<-^^ + 8VA|0| <-^<-ReA (13.33) 

2 A 

v'{u) = + 2n^u + 0{\n\u) + . (13.34) 

Since ReR is convex on (0, ttg] and has its minimum at uq, we know that the 
function ReR is monotone decreasing on the interval (0, uq]. This implies (113.32p . 
Moreover, using the results of the analysis in mm, we know that (113.18p holds with 
rigorous error bounds. This implies (113. 7|) for sufficiently large C. The inequality (113.61) 
follows from the estimate 
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(13.35) 


where we used (113.33p . (113.27p and (lll.2np . In order to prove (jl3.3ip . we again ap¬ 
ply (113.241) and estimate the two resulting integrals by 
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where in the last step we applied ()11.20p . (|13.33p and ()10.2I) . 
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^l(Ao) 



Figure 7. Annular regions. 


In order to estimate £, we proceed exactly as in the proof of Lemma [l3.7[ using the 
following estimates: 
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This concludes the proof. 


< 

rs-/ 


□ 


14. Annular Regions where |Cl| is Bounded Below 

We now construct regions in the complex A-plane in which |(^x,| is bounded below. 
For given Aq and AA > 0 we introduce the annular region 


Al{Xo) := {a gc| |Re(A- Ao)| G (AA,5AA) 
and |lm(A —Ao)| G (AA, 5AA)| 


(14.1) 


(see Figure ED. We again choose ui and Ur as the boundaries of the WKB regions 
(see (|lU. 8 p . (|1U.141) 1. With these choices, we have the following result. 

Proposition 14.1. For sufficiently large |n|, the following statement holds: Choose 
any e < 1/(1000). Suppose that for a given Aq, 

|a(l^max)| <e. (14.2) 

Then, choosing 

-Ur 


AA = 75e 






^0 


(14.3) 


it follows that 


e < |CL('Wmax)| < 1000e for all A G Al{\o) ■ 

The remainder of this section is devoted to the proof of this proposition. 
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14.1. Estimates in the WKB Region. We again consider the family of solutions 
of the Riccati equation (12.141) with initial conditions (|10.16p and (I11.13P (where uq 
is given by (111.121) or (|11.171) 1. We again let (p be the corresponding solution of the 
Sturm-Liouville equation normalized according to (16.71) . 


Proposition 14.2. Assume that 



(14.4) 


Moreover, assume that on [ui,Ur\ the WKB conditions in Proposition \ll.l\ as well as 
the inequalities (I13.6|) - (ll3.9p are satisfied. Then by choosing the constants Si,... S 3 
sufficiently large, we can arrange that the following statement holds. Choosing AA 
according to ()14.3p . for every A G Al{Xq), 


35 e < 




TOe < 


\U£ 

logK{X)\ll-log K{Xo)\ 


< 190 e 


Ui 


< 380e , 


with an arbitrarily small error. 


Proof. We analyze what the condition (114.4h means. Combining the estimate 
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with (114.3p . we conclude that 


AA < sup | 1 A| . 


^3 (ui,Ur) 


(14.5) 


This shows that varying A on the scale AA keeps the form of the potential V unchanged, 
up to an arbitrarily small error. 

By choosing the constants ci and C 2 sufficiently large, we can arrange that the as¬ 
sumptions of Lemma [13.4l and Proposition 113.51 hold. Possibly by further increasing ci 
and C 2 , we can arrange that on the interval [uQ,Ur], the WKB approximation 
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W- 


4K 


holds with an arbitrarily small error. Moreover, we can make the error term £ in 
Proposition 113.51 as well as the right side in (113.101) as small as we like. Hence on the 
interval {uo,Ur) and at A = Aq, 
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where ~ means “up to an arbitrarily small error.” 
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Since according to (|14.5p . the form of the potential is nearly constant on the 
scale AA, it follows that for any A G Ai(Ao), 


PI 

J un 


1 


luo 2lm\/V 
log A'(A)|;; - log A(Ao)|- O. - Im(A - A„) £ 

Using the form of the annular region, we obtain the result. 


-^0 




Re(A — Ao) 

32 AA 

^ Im(A - Aq) 
^ 16 AA 


□ 


14.2. Estimates of OxCl in the Airy and Parabolic Cylinder Regions. Before 
obtaining the estimates, we explain how we can arrange that the imaginary part of the 
potential satisfies the assumptions (a) or (b) in Proposition II 1.91 and Proposition lll.lTl 
To this end, we make use of the fact that taking the complex conjugate of the Sturm- 
Liouville equation (12.71) is again of Sturm-Liouville form, 

but now with the opposite sign of ImU. For the construction of the resolvent, we are 
free work either with the original equation or with the complex conjugate equation 
because if the resolvent of (jl4.6p has been constructed, the corresponding resolvent 
of (j2.7p is obtained simply by complex conjugation, preserving all our estimates. With 
this in mind, we can proceed as follows: In the parabolic cylinder region, the func¬ 
tion ImU either has a zero, or it is everywhere positive or negative. If it has a zero, 
we are in case (b) of Proposition 111.91 If it is everywhere positive, we are in case (a) 
of Proposition 111.91 If it is everywhere negative, we work with the complex conju¬ 
gate equation and are again in case (a). If we are in the Airy case. Lemma 111.121 
gives three possible cases. In case (i), we may apply Proposition 111.11] in case (a). In 
case (ii), we work with the complex conjugate equation and again apply case (a) in 
Proposition 111.111 Finally, in case (hi) the assumption (b) in Proposition 111.111 are 
satisfied for <5 < ^. We conclude that with this procedure, the assumptions (a) or (b) 
in Proposition 111.91 and Proposition (|ll.lip can always be satisfied. In what follows, 
we can take them for granted. 

We consider the interval [uj.,v\, where (cf. pO.lSp l 

{ u+ in the Airy case 

■Wmax in the parabolic cylinder case (14.7) 

Wmax in the WKB case . 

The following lemma is trivial in the WKB case because in this case the interval [ur,v\ 
reduces to a single point (cf. (|10.14p and (|14.7p ). But in the Airy and parabolic cylinder 
cases, the next lemma gives control of the region near the zero of Re V. 

Lemma 14.3. For any u G [ur,u], 

{dxCiiu) - dxCL{ur)\ 

< ^|<9Alog(/>(Ur)| {U -Ur) + Y |^Ay(Wr)| {u - Urf + (u - Urfj . 
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Proof. We first integrate the differential equation (17.121) from Ur to u to obtain 

J Uq" 

In the parabolic cylinder case, we now apply Proposition 111.91 Similarly, in the Airy 
case, we apply Proposition 111.111 This gives 

\(l>‘^{u)dxy{u)\ < \(l)‘^{Ur)dxy{Ur)\ + (u - Ur) (14.8) 

Next, we want to estimate dxlogf). Differentiating the relation dulogf) = y with 
respect to A, we obtain 


Integrating from Ur to u gives 

■ 

J Uj' 

Now we can apply (|14.8p and again Proposition 111.91 respectively Proposition 111.111 
to obtain 

|9Alog(^(u)| - I^AlogI < J 

< (^\(i)‘^iUr)dxy{Ur)\ (u - Ur) + y {u - Ur)^ 

g4 

= 02 |9A2/(ttr)| (n - n^) + Y (rt - . 

Finally, we compute the A-derivative of Cl, 

Using the above estimate for dxlog (j){u), we obtain 

9^2 pu 

\dxCL{u) - dxCliUr)] < I^Alog^l 

np2 ru / p4 \ 

< Y J _ (^|^Alog<?l(u,.)| +e| \dxy{Ur)\ (U - Ur) + Y (W - J ■ 
Carrying out the integral gives the result. □ 


AAiog0 = aA!/. 




Pit 


14.3. Proof of the Lower Bound for |^£,|. According to the estimates near the 
poles iLemmas 111.71 and 111. 8 p . by suitably increasing 65 we can arrange that 

\CLiue)\<S. (14.9) 

As is obvious from the expansion of the potential in (111.141) and (111.181) . by further 
increasing C 5 we can arrange that on [u£,no]. 


ReU>^>^|0|. 


2 
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Using furthermore that ImU is bounded by < |ri| (see again (|10.ip and ()9.3p l. pos¬ 
sibly by again increasing C 3 we can again arrange that the real part of the potential 
dominates the imaginary part. Hence 




ImU 

V\^v\ 


< {uo - Ul) 


1^1 < ^0 



(14.10) 


In view of the value of uq as given in (I11.12p and dll.lOj] . one concludes that by 
further increasing 65, we can make the left side in (|14.10p as small as we like. In view 
of the formula for the radius of the osculating circle in the WKB region (112.lOp . this 
means that that R is constant on the interval [ui,uq], up to an arbitrarily small error. 
Since R{uo) = ^ (as is again obvious from (112.101) 1. we can thus arrange by choosing 6 
sufficiently small that 


R{ue) 



(14.11) 


Combining (114.9|) with (jl4.1ip and using that Cl hes on the osculating circle with 
center p, the triangle inequality gives 




1 

2 


< 2e . 


Moreover, we know from (jl2.10l) that p is nearly constant in the WKB region. Thus 
we can arrange that 


IPI 




on [u£, Ur\ . 


Combining ()14.2p with the estimates of Lemma 111.131 we obtain the following esti¬ 
mates, 



Ci(^)L+ -Cl(Ao)|„^ 


3e 


> 

“Cl(Ao)|„^ 

- 


- 3e 



-b 

amii: - 

-b 3e 


In order to estimate |Cl(A)|^ — Cl(Ao)|^ , we go back to the osculating circle estimates 
of Proposition 114.21 Knowing that at ug, the function Cl is small (114.9p and the 
osculating circle has a radius close to one half (114.lip , the mean value theorem yields 


logKr(A)|„^ - logiL(Ao)U^ = logR{X)\ut -'^ogR{Xo)\ue 


R 


R{X)\^, - R{Xo)\ 


ue 


<,--e 


■2 5e 

2e<-, 


Moreover, elementary trigonometry shows that the angles satisfy the inequality 


1 


— e 2 sin 


m\u, - ^{Xo)\u, 


< 


|C(A)U,-C(Ao)| 


U( 


< 2s 


and thus 


^?(A)|n,-^?(Ao)U, 


< 2 arcsin 




-1 


e <5e. 
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In order to control the behavior on the interval [ui^Ur], we apply Proposition 114.21 
(note that the condition (I13.6P is satisfied in view of Proposition 112.1]) . It follows that 
one of the following two inequalities holds: 


30e<|'’P)L.-'’(V)Lj< 200E 
logA'(A)|„ -logif(A„)|„__ < 4006, 


60e < 

This in turn implies that the change of |Ci| can be estimated from below and above 
by 


Cl(A)|m^ - CL(Ao)k 
CL(^)|Mr “ CL(Ao)k 


> 10ei?(Ao)U, - p(A)|„^ - p(Ao)|n 

< 400e i?(Ao)|ur + P("^)kr ~ P(-^o)| 


> 4e 
< 500e . 


Finally, the estimate of Lemma 114.31 shows that the change of C,l in the Airy region 
or the parabolic cylinder region is much smaller than the change of Ql in the WKB 
region. In particular, by choosing |n| sufficiently large (and noting that the size of the 
parabolic cylinder and Airy regions tends to zero as |n| ^ oo), we can arrange that 


“Cl(Ao)| 

This concludes the proof of Proposition 114.11 


< e . 


15. The Green’s Function for a Double-Well Potential 

The goal of this section is to derive pointwise estimates of the Green’s function. For 
the statement of the result, we need the parameter u defined as follows. If Re V (umax) < 
0, we simply set u = Umax (this case includes the WBK case and part of the parabolic 
cylinder case). If conversely ReP(umax) > 0, we denote the zeros of the real part of 
the potential in the parabolic cylinder or Airy regions by 

Re V(vl) = 0 = Re R(vr) and vl < Umax < vn ■ 

Then u is defined by the requirement that 


ru rvji 

/ ReVv= / Re^R. 


(15.1) 


Since our potential is almost symmetric around the points u and Umax are all close 
to ^. Working with our definitions has the advantage that we do not need to quantify 
how close these points are. 

Here is the main result of this section. 

Proposition 15.1. Assume that for suitable constants C,6 > 0, the following inequal¬ 
ities hold, 

4 


|Cl(u) (/>i(u)|, \Cr{u) 4>ji{u)\ > 
\Cl\ <C on (0,u] , 


lyiiu) - yR{u)\ 

ICrI < G on [u,7r) 

|Cl(«)|, |Cr(«)| > 


Then, setting 


S = min 


(|</>L(«)|^ |</>r(«)|^) \ {yL-yR){u)\ 


(15.2) 

(15.3) 

(15.4) 

(15.5) 
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the kernel of the Green’s function sx{u, u') for 0 < m < u' < vr is bounded by 

1 


|sa(w,w')| ^ 


52 

C C 
+ 


5 

C C 
+ 


(5 52 s 


l^niu) 4>Riu')\ 


if u < u and u' > u 
if u, u' > 
if u, u' < 


1 /I 1 

+ T + 


1 


' V<5 ' 5‘^EJ \f>R\^ 

1 /I I \ r I 

+ T + 


We begin by estimating the Wronskian from below. 


\(t^Riu) ifu<u<u^ 

\(t)^{u)(l)L{u)\ ifu^<u'<u. 


Lemma 15.2. Assume that (jl5.2p holds. Then 


\wi 4 >LAR)\ > 77 Cl{u) Cr{u) w{(t)LAR) 


Proof. We evaluate (jl0.2ip at u = u, take the absolute value and use the inequali¬ 
ties (I15.2p . This gives 

^ </>i?CR {VL - yR)\iu) = ^ \CLCRwi4>L,4>R)\iu) , 

concluding the proof. □ 


In order to estimate the Green’s function, we need to control both 0° and on 
the whole interval (0,7r). The estimates so far, however, only give us control of on 
the interval (0, (and similarly of (p^ on the interval [f ,vr]). The following lemma 
gives a formula for (jf^ on the remaining interval [f ,vr]. 


Lemma 15.3. For any tt G [^,7r], 


4>^{u)-(t>R{u){cL{u) +Cr„ hm V 

V (Pr{u) u v\ 0 (j)L[v) J 


(15.6) 

where 



y 4 >r{v) (pRiu) ,, , \r 

A ! \ = A i'-\ +'>^WL,(Pr)Cl{u) . 

9 L[V) (Pl(u) 


(15.7) 

Proof. First, 



= +Mn) r {*})'( 

(Pr{u) 9 r[u) Ju ^(pR^ 

t) dr . 



Computing the derivative on the right gives 

^ ^ w{4)2,(Pr) 

Integrating this equation from u to u, we obtain 

4>Liu) = 4>Riu) + wifAl.cpR) (pniu) f ^ . 

(pRiu) Ju <pR 
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The Wronskian appearing here is most conveniently computed asymptotically at the 
origin, 

w{(I)lAr) = lim {{(I>^)'{v)^r{v) - {(t>R)'{v)) = lim , 

v\0 v\0 (pL\V) 

where in the last step we differentiated (j2.17j) and used the asymptotics of the funda¬ 
mental solutions as stated in (j2.11j) . (j2.12j) and (I2.13h . Applying (j6.2h gives (jl5.6j) . 

In order to prove (jl5.7ji . we begin with the computation 


j. ! \ j, 1 ^ j. ! ^ ^R^R 

(t>R[v) = (Pl{v) , . = 4>l{v) 


= 4>Liv) 


(Piiv) 

(t>R{u) 


Mu) 


^4>l 


+ 


Mv) f 

J V 


wi4>L,4>R) 


(j)L{u) ' Jv 4>i 

We now divide by (I)l{v) and take the limit v \0. This concludes the proof. 


□ 


Proof of Provosition \15.1\ Applying Lemma 115.21 our task is to estimate the absolute 
value of the expression 


E := 


(pLiR)(pRiu' 


Cl{u)Cr{u)w{4)lAr) 


for all 0 < u < < TT 


(15.8) 


We consider the different cases after each other. In the case u < u and u' > u, we 
apply (115.41) to (jl5.8p to obtain 


IPI < J_ 

“ (52 w{4>lAr) 


Using that 

\w{(t>LAR)\ = \(t)L{u)(t)R{u) {VL - yR){u)\ > S , 
we obtain the desired estimate. 

Using the symmetry under reflections at it remains to consider the case u, u' > u. 
Applying Lemma ll5.3l a straightforward computation yields 


T 7 ? Cr(,R^ Cr(.R) ± \ j, ( l\ /U f l\ 

E = --WTT- 4>r{u) f)R{u ) Cr{u ) 


CRiu) 


1 


CRiu) - CRiu) 


(pRiu) (t>Riu') CRiu') 


VLiu) - VRiu) CLiu)(j)Liu)‘^ CRiu) 

— 7m - 7m /- fMA (M2 C^Ri'^) C>Riu') CRiu') . 

VLiu) - VRiu) Cr{u) (pRiu)^ 


(15.9) 

(15.10) 

(15.11) 


In the case u > u^, we estimate these terms by 


IdlMllI < J \(j)Riu) (f)^iu')\ 
|(ll5Jilj)| < ^ \(j)Riu)(f^iu')\ 
|(ll5JTj)| < ^ \(j)R{u)(fRiu)\ 


giving the desired estimate. 
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In the remaining case u < u < u!^, we use the identity 

ru ^ 


to obtain 


ICniu) - Cr(^^)| = 


lu 4 >r 


1(115^1 < T |^r(^)</>rK)| 


In (t>R 


Ips.iupl < ^ |<^R(n) (j)R{u')\ 

lOEIDI < |'#'R(«)<^gM| ■ 


III 4>r 


This gives the result. 


□ 


16. Deforming the Potential 


In the following estimates, we distinguish the cases that the real part of the po¬ 
tential is positive and large near u = Umax or that it is negative or small there. We 
refer to these cases as the double-well case and the single-well case, respectively. 
Qualitatively speaking, in the double-well case the potential looks like in the Airy case, 
whereas the single-well case comprises the WKB and parabolic cylinder cases. How¬ 
ever, these regions are not exactly the same, making it necessary to use a new notation. 
The important difference is that the single- and double-well cases are defined without 
referring to the constants Ci, 64 ,.... Instead, we work with the integrals (jl5.ip and 
introduce a new constant X > 0 (which will be specified Section [16. ll below). More pre¬ 
cisely, if Re V (umax) = 0, we are by definition in the single-well case. If Re V (umax) > 0 
(so that u is defined by (115.111 b we distinguish the 


single-well case : 

< 

double-well case : 



ReVv<X 

ReW>X. 


(16.1) 


16.1. Estimates in the Double-Well Case. 


Lemma 16.1. By choosing X sufficiently large, we can arrange that in the double-well 
case the following estimates hold: 


\(j)L{u)\, \(Pr{u)\ > 

(16.2) 

\yLiu) - yRiu)\ > y/\V{u)\ . 

(16.3) 

Moreover, the function S defined bv (115.51) is bounded bv 


S > exp ^2 J Re . 

(16.4) 


Proof. If X is chosen sufficiently large, we know from the estimates in Section 111.61 
that the function (and similarly \4>r\) can be approximated by the WKB wave 
function ()11.39p . Moreover, as stated after (lll.39p . the coefficient of the exponen¬ 
tially increasing fundamental solution is non-zero. This implies ()16.2p . Moreover, 
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differentiating (111.391) . one finds that yL{u) ~ ^/V(u). Similarly, yR{u) ~ --s/V(u), 
proving (jl6.3l) . 

Finally, the estimate (jl6.4l) follows immediately by using (116.2p and ()16.3I) in (|15.5p . 

□ 


From now on, we choose X so large that the the statement of this lemma applies. 


Lemma 16.2. For any e > 0, by increasing % we can arrange that for every eigen¬ 
value Ao, 

|Cl(w)|<| or |Cr(«)|<|- (16.5) 

Proof. Let Aq be an eigenvalue. Then the eigenvalue condition (I10.22P is satisfied 
at u = u, implying that 

1 1 , , . 

+ 1,2 A I ^ IVL - VRl at u . 

\wlCl\ IwrCr] 

Applying Lemma ll6.ll we obtain 


1 ^ 1 
lawL^ KMOl 



> e 


23C 


Applying again the argument after (jll.lOp . by increasing Ci we can make the ex¬ 
ponential factors as small as we like. Hence, for any given e > 0 we can arrange 
that 

1 1 2 
|Cl(w)I \CRiu)\ ~ e ■ 

This gives the claim. □ 


We now apply Proposition 114.11 In order to combine estimates for (r and Cr, we 
modify (114.3p according to 


A A := 75 e min 





(16.6) 


Thus we choose AA so small that we can work with the same AA both for Qr and Qr. 
Decreasing AA in this way can be described equivalently by making the parameter e 
in Proposition (jl4.ip smaller. The reason why this procedure is unproblematic is that 
the parameter e is changed at most by a uniform constant: 


Lemma 16.3. 



Proof. Similar as in the proof of Lemma 112.121 and Section 112.51 we make use of the 
fact that the potential is approximately reflection symmetric around u = ^. We thus 
obtain 





|H(u) — V{7r — u)\ 

ml 


1^1 

CilDI 



1 


This gives the result. 


□ 
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^l(Ao) 



Al(Ao) 




\ 

\ 

^ — 







Figure 8. Choice of the closed contour F. 


Lemma 16.4. Let Aq he an eigenvalue. Then, under the assumptions of the above 
Lemmas \16.1\ and \16.3l. there is a closed contour F of length at most 20 AA such that 

|CLUCR|>e along F. 

Proof. Let Aq be an eigenvalue. Using Lemma 116.21 we can arrange that (|16.5I) holds. 
Without loss of generality we may assume that 

l&(«)| < I, 

because otherwise we repeat the proof with L and R interchanged. Applying Propo¬ 
sition [TTTl we know that |CL(fi)| > £ inside the annular region Al{\q). We choose a 
contour F inside Al{\q) close to the inner boundary (see the left of Figure [8|). If the 
inequality |CR(fi)| > ^ holds along the contour F, there is nothing left to prove. Other¬ 
wise, there is a point u on F such that |CR(fi)| < Again applying Proposition 114.11 
it follows that |Cr('w) > £ inside the annular region A/j(u). This makes it possible to 
choose the contour F as shown on the right of Figure [H □ 


Proposition 16.5. The resolvent sx exists along F and is bounded by 

e |(/)°(u) (()^(n')| if u < u and u'> u 


if u, u' > 
if u, u' < 


|(/)R(u)(/)g(n')| ifu<u<u^ 

. |(()°(n)(()L(ri')| if < u'< u . 


Proof. We apply Proposition II5.11 By symmetry, in the estimates it suffices to consider 
the index L. 

We begin by estimating |^£,|. In the WKB region, the total variation of the ra¬ 
dius R of the osculating circle is under uniform control by Proposition 112.11 In the 
region (0, ui) near the pole, on the other hand, this total variation is even small by 
Lemma 112.61 In the Airy region. Proposition 111.111 provides uniform estimates. Fi¬ 
nally, Lemma [T lT 3] controls the total variation of Ql on the interval [n+, u]. Combining 
these results, we conclude that 


IClI < 1 on(0,ii]. 


This gives ()15.3p . Moreover, we know from Lemma 116.41 that ()15.4I) holds for 6 = e. 
We conclude that the assumptions ()15.3I) and (|15.4I1 in Proposition 115.11 are satisfied. 
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Let us verify the assumption (I15.2I1 . From (|16.2p and (I16.3P we know that 

^2X 

y/\V\ V JvL ' y ^ IVLiu) - yRiu)\ ■ 


\Cl{u)(PI{u)\ > -^= exp (2 J ReW 


Thus by further increasing % we can arrange (I15.2p . Therefore, Proposition 115.11 
applies. 

We estimate the integrals over l/*/*^ by 


= - rAe-2/5Rex/y 


Ju 

< ( p-2/“^ Re^W _ 2/;?Re^W^ < -2£«Re^W 




Using the lower bound for S in (jl6.4p as well as (jl5.ip gives the desired estimate. □ 

16.2. Estimates in the Single-Well Case. In the previous section we chose the 
parameter JC. We now proceed with estimates in the single-well case for this fixed 
value of JC. 

Proposition 16.6. There is a suitable choice of e < 1/(1000) such that for sufficiently 
large Ci the following statement holds. If Aq is eigenvalue in the single-well case, then 
the resolvent exists inside the annulus Al{Xq) (defined by (114. ip ) and is bounded by 

\sxiu,u')\ <c\(j)Y.iu)4>Riu')\ ifu<u' 

(where the constant c may depend on JC and s). 

Proof. Let Aq be an eigenvalue. Then the eigenvalue condition (|10.22p holds for any u E 
(0, tt). We now choose a specific value u where the analysis of the expression in (jlO.221) 
is particularly simple. To this end, we consider the osculating circle corresponding 
to 4>l on the interval (0, Ur). According to Lemma fl2.6l . we can make |C(i^^)| arbitrarily 
small. Thus for any 5 > 0 by increasing Ci we can arrange that 

\ R { ui ) - |p(rt£)|| < d . 

Next, in the WKB region {ug,Ur), we know from Proposition 112.71 that the center p 
of the osculating circle is approximately fixed, whereas the change of the radius R is 
given explicitly in terms of the integral of Re\/U. Using the estimate (jl2.21|] together 
with Proposition 112.21 for any given <5 > 0 we can arrange that that 

\R{ur) - |p(nr)|| < 6 . 

If we start at u = Ur and decrease u, the function C(rt) will move along the osculating 
circle with an angular velocity as given by (16.9p . We choose u as the largest value of u 
where the angle d is such that ^ is close to the origin, i.e. 

|a(«)|<25. (16.7) 

In what follows, we evaluate all functions at the point u. For ease in notation, we 
shall omit the arguments of these functions. From the estimates in Section [TT] we know 
that 




1 


vT^ 


(16.8) 


(where the constants in the upper and lower bounds may depend on IK). Moreover, 

IVL - VrI < VWl ■ (16.9) 
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Therefore, the eigenvalue condition (I10.22P gives rise to the estimate 


ICrI 


1 

IClI 


\yL - yR\ 


> 


implying that 


|Cr| <<5. 


(16.10) 


Combining (|16.7p and (116.lOp we conclude that for any 5 > 0, we can arrange by 
increasing Ci that at u, 

\Cl\, |Cr| < ^ • 

We now vary the spectral parameter A (for fixed u and 3C). Our goal is to analyze 
how the left side of ()10.22p depends on A. To this end, it is most convenient to include 
the phase of the factor ((P‘ into the corresponding factor Thus, using the notation 
in (I63|) , we set 

^ gQ ^ 2 ^ _ |^|2 ^ ^ 

Then, according to (j6.6p . 

Hence (I10.2ip can be written as 




K 








= yL-yR + 


It follows that 


w{4P 


L> 


\(t>R?ZR 

\(I>R?Zr - \4)l?zl 


= yL-yR + 


\4>R?Zr - \(t>L?ZL 

\(t>L\^ZL \(t)R\^ZR 




> 


\yL - yR\ 


\(t>L?ZL \(t)R\^ZR 

Now we can vary the angle 'd and the radius R = of the osculating circles using 

the formulas in Proposition 113.51 Keeping in mind that the variation of 'dR and Kr 
involves a minus sign (because we consider the differential equations backwards in u), 
one sees that the variations of Re zr and Re zr (and similarly of Im zr and Im zr) have 
opposite signs. 

We now choose A in the annulus Ar{Xq) as dehned by ()14.ip . The resulting vari¬ 
ations of zr and zr are obtained by integrating the infinitesimal variations, exactly 
as explained in the proof of Proposition 114.11 Given e > 0, by choosing 5 sufficiently 
small we can arrange the variation is much larger than 2 ; at A = Aq (see the formulas of 
Proposition 114.21) . Thus the above consideration for the sign of inhnitesimal variations 
implies that RezR and Rezi? (and similarly \m.ZR and Imz/j) have opposite signs. 
This gives rise to the inequality 

\\4>R?zr - \(I)r\‘^zr\ > \(I)\zr\ . 

It follows that inside the annulus Ar{Xq), 


w{(I)r^(Pr) 




> 


\^rCr\ 


\yL - yR\ 


Using ()16.8I) and (116.9p together with the fact that |Cl| and \(r\ are uniformly bounded 
from above, we conclude after choosing e sufficiently small, the first summand ma¬ 
jorizes the second summand, i.e. 


w{4'r^4> 




> 


1 


\^rCr\ ' 
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Multiplying this inequality by and we obtain 




\4>r\ 


Using again that = \zl\ is bounded from below, we conclude that 
Finally, we use this estimate in (j2.19p . This concludes the proof. 


□ 


16.3. Tracking the Eigenvalues. We are now in the position to complete the proof 
of Theorem 0 In order to locate the eigenvalues, we use the following deformation 
argument. We consider for r € [0,1] the homotopy 

Vr = ReU + rlmU 

with V as in (12.7p . Then at r = 0, the potential is real. As a consequence, the 
Sturm-Liouville operator is self-adjoint. It has a purely discrete spectrum with real 
eigenvalues. If r > 0, on the other hand, the potential is complex. Our method is to 
track each eigenvalue as r is increased. 

At T = 0, we choose the contours as shown in Figure [5l with N chosen as follows. 
First, we choose N at least as large as in Proposition 17.71 for C 4 so large that the 
eigenvalue Ajv satisfies the inequality (jl0.2p . The next lemma makes it possible to 
arrange a spectral gap. 

Lemma 16.7. By increasing N at most by four, we can arrange that for a suitable 
constant c > 0 and large |n|, 

Atv — Atv —1 ^ c|n| . (16.11) 


Proof. In preparation, we want to show that there is a constant Ci > 0 such that for 
sufficiently large |n|, the eigenvalues of the Hamiltonian (17.2p satisfy the inequality 


Xn+4 — Aat > Cl |n| . 


(16.12) 


To this end, we first note that the analysis in the proof of Proposition 17.71 shows that 
the solution has no zeros on the interval {0,u^) near the pole at tt = 0. Similarly, 
the solution (j)^ has no zeros on the interval (u^, tt) near the pole ai u = ir. Moreover, 
on the interval where the potential is non-negative, we know from Lemma [7.4l 

that (fP has at most one zero. Therefore, counting the number of zeros of cj)^ on the 
intervals and this number differs at Xn +4 and Xn at least by three. 

As a consequence, on one of the intervals {u^,u^) or {u^,u^), the number of zeros 
of (fP differs at least by two. By symmetry, we may assume without loss of generality 
that this is the case on the interval (ti^,ni). It follows from Lemma 17.51 that 


(^Imy|Ajv+4-ImyUjv) 


> vr . 


Applying the mean-value theorem in the parameter A, we obtain 


/■“+ 

{Xn+4 - Xn) sup / \dxy\>7r 

AG[Ajv,Aiv_l_4] Ju^ 
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In order to estimate d\y, we use the formula (|7.13p . Inserting the asymptotic expan¬ 
sions near the poles (see (I11.15P and (|11.221) 1 and using the rigorous estimates in [121 
Section 5 and 8], a straightforward computation gives the inequality (116.121) . 

The inequality (|16.12l) shows that for any 11, one can choose N such that (116. lip 
holds (for c = Ci/4). In order to show that N can be chosen uniformly in 11 (for 
large |11|), we make use of the fact that the leading powers in H in the potential (12.811 are 
symmetric around u = ^. This implies the small eigenvalues have the same asymptotics 
as ReH — )• ±oo (with the corresponding eigenfunctions transforming as 4>n{u) —>■ 
(^n(7r — n)). As a consequence, we can choose N such that (116.12 p holds for all H 
with |11| sufficiently large. □ 


For the chosen iV, we choose the contour Tg such that it encloses the lowest N 
eigenvalues, where N is chosen at least as large as in Proposition 17.71 The other 
contours enclose one eigenvalue (see Figure [5]). Then all eigenvalues not enclosed 
by Tg satisfy the inequality (110.21) . 

Now we continuously change the parameter r and follow each of the eigenvalues. 
We also continuously deform the contours Fg, Fi, ... such that they enclose the cor¬ 
responding points in the spectrum. If Ag is an eigenvalue in the single-well case, 
according to Proposition 116.61 the resolvent is well-defined and bounded along a closed 
contour P in the annular region Aj;,(Ag) (which can be chosen for example as on the left 
of Figure [8]). For an eigenvalue Ag in the double-well case, on the other hand, there are 
two possible subcases. Either the eigenvalue is isolated in the sense that we can again 
choose the contour as on the left of Figure [71 Or else the eigenvalue can be close to 
another point in the spectrum, in which case we choose the contour as on the right of 
Figure [HI In both subcases, the resolvent estimate of Proposition 116.51 applies. Finally, 
we enclose the remaining N lowest spectral points again by a contour Fg. Dehning the 
corresponding operators Qn by 


Qn 



neNg, 


we obtain a family of operators. 

For clarity, we point out that spectral points may move from the single-well case to 
the double-well case and vice versa. Moreover, the contours F^ do not need to be chosen 
continuously in r. Indeed, if two eigenvalues come close together, the corresponding 
contours must be changed discontinuously because the two contours enclosing the two 
eigenvalues (as on the left of Figure |8|) must be joined to form one contour (as on the 
right of Figure Ej). As a consequence, the operators Qn will in general not depend 
continuously in r. 

It remains to verify a-posteriori that the inequality ()in.2p holds for all r € [0,1] 
and for all spectral points not enclosed by Fg. To this end, we first point out that 
the eigenvalues may change considerably compared to the size of the gaps between 
neighboring eigenvalues (this is why the theory of slightly self-adjoint perturbations 
does not apply here). But the previous methods tell us how the eigenvalues change 
in T (in particular see Section (U Section [10.21 and Section [12.3p . In simple terms, 
these results show that the eigenvalues must satisfy the complex Bohr-Sommerfeld 
condition (|3.6p with well-defined errors. Combining this formula with the mean value 
theorem 


< 


IlmFl 
sup — - - 

re[0,l] V |lr| 




(16.13) 
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and using that on the interval {ui,Ur) the absolute value of the potential is larger 
than Cl |f7|, whereas | Imy| < |n|, one sees that the eigenvalues change in such a way 
that (jl0.2l) remains valid. 


16.4. Uniform Boundedness of the and Completeness. It remains to verify 
that the constructed operators Qo,Qi,... at r = 1 have all the required properties. 
The idempotence and mutual orthogonality of these operators follows immediately 
from Lemma E21 We now proceed by estimating the sup-norm of the operators Qn- 
For the operator Qq we again apply the theory of slightly non-selfadjoint perturbations: 

Proposition 16.8. There is a constant C 2 such that for all sufficiently large |n|, 

IIQoll < C2 . 

Proof. In view of the gap estimate (|16.11l] . by increasing |fl| we can arrange that | Im V\ 
is much smaller than the gap. This makes it possible to find a contour Tq enclosing 
the first N spectral point whose distance to the spectrum is much larger than | ImU|. 
This makes it possible to estimate the corresponding contour integral in (18.4p for n = 0 
by estimating the Neumann series (18.31) . This gives the result. □ 

We next estimate the operators Qi with i > 1. We begin with a preparatory lemma. 


Lemma 16.9. Along the contours Ti, r 2 ,..., the Hilbert-Schmidt norm of the resolvent 
is bounded by 


I'SaIIhs ^ 


r- 1 

1 

luf 

o 

4 


-^0 


(with Aq as in (116.6|) ). 


|sa||hs — 2 
= 2 


Proof. The Hilbert-Schmidt norm of the resolvent can be expressed in terms of its 
kernel by 

P7V P7V 

\\s\\\m= du du'\sxiu,u')\‘^ . 

Jo Jo 

We begin with the single-well case. Proposition 116.61 gives the estimate 

f du f du'\(P l{u) 

Jo J u 

f du f du'\CLiu) (piiu) Cr{u')(P^{u')\'^ . (16.14) 

Jo J u 

Near the pole at = 0, the functions cfi and 4)r may have a pole (see ()11.22p l. On 
the other hand, the function Cl{u) vanishes at rt = 0. As a result, the integrand 
in (116.141) is bounded near the poles, as the following argument shows. We introduce 
the functions pL and pR by 

Pl{u) = sup IClI and pr{u) = sup \Cr\ . 

Then the integrand in (|16.14l) can be bounded by 

|Cl(^^) 4>l{u) Cr{u')4>^{u')\‘^ < pl{u) \4>l{u)\'^Pr{u) pl{u)\( j)R{u')\^ pr{u) , 
giving rise to the estimate 

PIT 

PaIIL^ n / Pl{u)\(I)c{u)\^ Pr{u) . 

c=L,Rd(> 


(16.15) 
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In order to estimate the obtained integrand near u = 0, we first note that Cr is 
uniformly bounded. Considering the asymptotics near u = 0 (see (|11.16p and (|11.22l) l. 
one sees that the function \4>c\‘^ may have a pole at u = 0. But in this case, inserting 
the asymptotics into (|6.ip . one finds that |Ci,| (and therefore also pi) tends to zero at 
the inverse rate. We thus conclude that the integrand in (116.15p is indeed bounded 
near the poles. 

Away from the poles, we can use the asymptotics (|16.8I) to obtain the result. This 
concludes the proof in the single-well case. 

In the double-well case, we work similarly with the resolvent estimate of Proposi¬ 
tion [T63J The behavior near the poles is estimated just as in the single-well case. The 
only additional issue is the behavior on the interval [u^,u^] in the Airy case. In this 
case, combining the exponential factor exp(—2 Re W) with the WKB asymptotics 
gives rise to the estimate 


e \<Pr{u)^r{u’)\ 




,-2/, 


minlti ,v 


Re\/T 


A straightforward computation using the exponentially decaying factor on the right 
gives the result. □ 


Proposition 16.10. There is a constant C 2 such that the operators Qi,Q 2 -, ■ ■ ■ «Fe 
hounded by 

11 Qn 11 ^ C2 , 


uniformly in 11. 


Proof. We estimate the contour integrals by 

IIQnll < A(r„) sup ||sa|| , 

AGFtt, 

where L{Tn) denotes the length of the contour. The length of the contour is bounded 
by 20 AA with AA as given by (|16.6I) (this is obvious in the single-well case, whereas 
in the double-well case it was proven in Lemma 116.41) . Applying Lemma 116.91 and 
Lemma 116.31 and using that the Hilbert-Schmidt norm majorizes the sup-norm, the 
result follows. □ 


It remains to prove completeness in the sense of (II.7|) with strong convergence of 
the series. Since completeness is a statement for fixed 12, we can rely on the theory of 
slightly non-selfadjoint perturbations. Namely, Proposition 18.11 yields that there is N 
(which might be much larger that the parameter N in the statement of Theorem II.Ih 
a family of operators Qo,Qi, ■ ■ ■ with the completeness property 

CO 

Qfi = 1 with strong convergence . 

n=0 

The spectral decompositions (Qn) and {Qn) are related to each other simply by forming 
finite sums, i.e. 

E<?» = E Qn (with A, A c No) 

ngA 

whenever the operators on the left and right describe the same spectral points. In 
particular, for large n, the operators Qn are sums of one or two of the operators Qn- 
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This implies that the series in (ll.7p also converges strongly. Since every the spectral 
point is taken into account in exactly one of the operators Qn, it follows that 

OO OO 

Ee» = E Qfi with strong convergence . 

n=0 n=0 

This concludes the proof of Theorem 11.11 
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